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ABSTRACT 

From high-resolution images of 23 Seyfert-1 galaxies at z=0.36 and z=0.57 obtained with 
the Near Infrared Camera and Multi-Object Spectrometer on board the Hubble Space Telescope 
(HST), we determine host-galaxy morphology, nuclear luminosity total host-galaxy luminosity 
and spheroid luminosity. Keck spectroscopy is used to estimate black hole mass (Mbh)- We 
study the cosmic evolution of the MBH-spheroid luminosity (L S ph) relation. In combination with 
our previous work, totaling 40 Seyfert-1 galaxies, the covered range in BH mass is substantially 
increased, allowing us to determine for the first time intrinsic scatter and correct evolutionary 
trends for selection effects. We re-analyze archival HST images of 19 local reverberation-mapped 
active galaxies to match the procedure adopted at intermediate redshift. Correcting spheroid 
luminosity for passive luminosity evolution and taking into account selection effects, we determine 
that at fixed present-day V-band spheroid luminosity, Mbh/ L^oc (1 + z) 2 8±1 - 2 . When including 
a sample of 44 quasars out to z = 4.5 taken from the literature, with luminosity and BH mass 
corrected to a self-consistent calibration, we extend the BH mass range to over two orders of 
magnitude, resulting in Men/^sphoc (1 + z) 1 °' 2 . The intrinsic scatter of the relation, assumed 
constant with redshift, is 0.3±0.1 dex (<0.6 dex at 95% CL). The evolutionary trend suggests 
that BH growth precedes spheroid assembly. Interestingly, the M^H-total host-galaxy luminosity 
relation is apparently non-evolving. It hints at either a more fundamental relation or that the 
spheroid grows by a redistribution of stars. However, the high-z sample does not follow this 
relation, indicating that major mergers may play the dominant role in growing spheroids above 
z ~ 1. 
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1. INTRODUCTION 

Supermassive black holes (BHs) seem to be 
ubiquitous in the center of spheroids - ellipti- 
cal galaxies and classic a l bulg e s of spirals (e.g.. 



Kor mendv fc Richstond Il995t iFerrarese fc Ford 



2005h . In the local Universe, tight empirical rela 



tions have been found between the mass of the BH 
(Mbh) and the properties of t he spheroid, i.e. stel- 
lar v elocity dispersion a (IFerrarese &: Merritt 



2000; Gebhardt 



et al 



2000). stellar mass (e.g., 
and luminosity (e.g 



Mar coni fc Hunt] l2003h . ,„ 

Haring fc Rixll2004h . The tightness of these rela- 



ogy, Stanford, CA 94305; rdb@slac.stanford.edu 
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tions is surprising, given the very different scales 
involved - from accretion onto the BH (//pc 
scale), the dynamical sphere of influence of the 
BH (pc scale) to the size of the spheroid (kpc 
scale) - and poses a challenge to any theoreti- 
cal model explaining their origin. In general, the 
correlations are believed to indicate a close con- 
nection between galaxy formation and evolution 
and the growth of the BH. A variety of theo- 
retical models have been developed to explain 
the observed relations, involving galaxy merg- 
ers and nuclear fe edback through quenching of 
star formation (e.g., 



these problems in distinct ways, e.g. by using the 
[O Hll emission lin e width as surrogate of a (e.g., 
Shield s et al 1 l2003h . or by using gravitational lens- 



ing t o super-resolve th e host galaxies of quasars 



2004; 



Peng et al 



Woo et al 



I2006alibh. Our group (llreu et al 
2006HTreu et al.ll2007tlwoo et al 



Volontcri ct al. 2003; 



Cauffmann fc Haehneltl feoOO: 



Ciotti fc Ostrikcr 2007: Hopkins et 



20071: iDi Matteo et al.l2008HHopkins et al.l2009bl) . 



Measuring the evolution with redshift of these 
correlations constrains theoretical interpretations 
and provi des important insights into the i r ori- 



gin (e.g.. ICr oton 



20061 : iRobertson et all 120061 



Hopkins et al.ll2007T ). For quiescent galaxies, the 
biggest challenge is to measure the BH mass, given 
the pc-scale sphere of influence of the BH which 
needs to be resolved sp atially through either gas 
or stellar dyn amics (see iGultekin et al.l 120091 and 



Graham 2008 for a recent com pilation and refer- 



ences therein; for a review see iFerrarese fc Ford 



2005 and references therein) or from X-ray spec- 



troscopy probing the existence of a central tem- 
perature peak of the inter s tellar medium (e.g., 
Brighenti & Mathews 1999c iHumphrev fc Buote 
2006t iHumphrev et all 120081 ). With current tech- 



nology, direct quiescent black hole mass measure- 
ments are thus limited to nearby galaxies. 

For active galaxies, for which nuclear luminos- 
ity is comparable to or larger than that of the 
host galaxy, the situation is virtually the oppo- 
site. Estimating BH masses within a factor of 
2-3 is fairly straightforward through empirically 
calibrated relations based on spectroscopic data 
measuring the kinematics of the broad-line re- 
gion ( BLR) (e.g.JWandel et al.lll999tlWoo fc Urrv 



2002; IVestergaardl l2002t IVestergaard fc Peterson 
2006t iMcGill et all 120081 ) . Unfortunately, the ac- 



tive galactic nuclei (AGN) often outshines the 
host galaxy, making it difficult to disentangle 
nuclear and host-galaxy light for an accurate 
measurement of the spheroid luminosity. Also, 
measuring a from stellar absorption lines is ham- 
pered by the contaminating AGN continuum and 
emission lines. Different groups have tackled 



2008) has focused on Seyfert-1 galaxies - for which 
the nucleus is not as bright as for quasars - at 
moderate redshifts (z — 0.36 and z = 0.57, cor- 
responding to look-back times of —4-6 Gyrs). 
The non-negligible stellar light produces strong 
enough absorption lines to m easure a from un- 
reso lved spectr a as shown by llreu et al. ( 2004 ) 
aflnd IWoo et al.l (|2006L 120081 . hereafter Paper I & 
III). At the same time, high resolution Hubble 
Space Telescope (HST) imaging allows for an ac- 
curate determination of the AGN luminosity (for 
an unbiased estimate of nuclear luminosity and 
hence Mbh) and sphe roid luminosity (to create 
the MeH-^sph relation; Treu et al. 20071 . hereafter 
Paper II). We are thus able to simultaneously 
study both the Mbh-c and MeH-^sph relations, 
allowing us to distinguish mechanisms causing 
evolution in a (e.g., dissipational merger events) 
and I/ sp h ( e -g- through passive evolution due to 
aging of the stellar population, or dissipationlcss 
mergers). 

Results presented in Paper I, II, and III sug- 
gest an offset with respect to the local relation- 
ships, which cannot be accounted for by known 
systematic uncertainties. At a given A/bh: in 
the range 10 8 -10 9 Mq, spheroids had smaller ve- 
locity dispersion and spheroid mass 6 Gyrs ago 
(z ~ 0.57), consistent with recent growth and evo- 
lution of intermediate-mass spheroids. Paper II 
concludes that the distant spheroids have to grow 
by -60% in stellar mass (AlogM sph = 0.20 ± 
0.14) at fixed black hole mass in the next 4 billion 
years to obey the local scaling relations if no signif- 
icant BH growth is assumed, consistent with the 
relatively low Eddington ratios. Indeed, the HST 
images reveal a large fraction of merging or inter- 
acting systems, suggesting that gas rich mergers 
will be responsible for the spheroid growth. 

Although tantalizing, the results presented in 
our previous papers suffer from several limitations. 
Samples were small, and the local comparison 
sample of Seyferts measured in a self-consistent 
manner was even smaller than the distant sam- 
ple, thus contributing substantially to the over- 
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all error budget. The limited range in black hole 
mass was insufficient to determine independently 
the offset of the scaling relation and its scatter, 
while taking into account selection effects. If the 
Mbh-u and MBH-L sp h relations of active galax- 
ies were not as tight as those for quiescent ones, 
selectio n effects could be mimicking evolut i onary 
trend s (|Treu et al.ll2007t lLauer et alJl2007t IPengl 
2007h . 

To overcome these limitations, we have now 
doubled the sample size (from 20 in Paper II to 40 
total here) and expanded the covered range of BH 
masses to lower masses (from log Mbh/Mq — 8 — 
8.8 in Paper II to log M B n/M Q = 7.5 - 8.8 here). 
We focus on the resulting BH mass - spheroid 
luminosity relation. The BH mass - a relation 
will be presented in a separate paper (Woo et al. 
2009, in preparation). We also analyze archival 
HST images of the sample of local Seyferts with 
reverberation-mapped (RM) A/bh in the same way 
as our intermediate-z objects, to eliminate possi- 
ble systematic offsets. Finally, we combine our 
results with data compiled from the literature and 
treated in a self-consistent manner to extend the 
redshift range over which we study evolution. For 
conciseness the three samples will be referred to 
as "intermediate-redshift" sample, "local" sample, 
and "high-redshift" sample, respectively. 

The paper is organized as follows. We sum- 
marize the properties of our intermediate-redshift 
Seyfert sample, observations, data reduction, and 
analysis in § [2j [3j and |4j § [5] summarizes the 
derived quantities, including the derivation of 
Mbh from Keck spectra. In §[6j we describe the lo- 
cal comparison sample consisting of reverberation- 
mapped AGNs, re-analyzed here, as well as the 
high-redshift comparison sample taken from the 
literature, calibrated for consistency with the 
other samples. We present our results in § [3 in- 
cluding host-galaxy morphology and merger rates, 
the evolution of the MeH-^sph relation, a full dis- 
cussion and treatment of selection effects, and a 
relation between BH mass and host-galaxy lu- 
minosity. We discuss the possible implications of 
our findings for the origin and evolution of the BH 
mass scaling relation in §[8] A summary is given in 
§0 In Appendix [A] we describe Monte Carlo sim- 
ulations used to probe our analysis and determine 
errors. Appendix [B] discusses the choice of the 
Sersic index in the adapted 2D surface-brightness 



fitting procedure. Details on the re-analysis of the 
HST images of the local RM AGNs are given in 
Appendix [C] Throughout the paper, we assume 
a Hubble constant of Hq = 70kms _1 Mpc^ 1 , J7a 
= 0.7 and I^ m = 0-3. M agnitudes are given in the 



AB system (|Okelll974l) 



2. SAMPLE SELECTION 

The selection of the sample of intermediate- 
redshift Seyfert- 1 galaxies is similar to the one in 
Paper I, II, and III, with the goal to extend the 
sample to (a) lower BH masses at z ~ 0.36 and (b) 
higher redshifts of z ~ 0.57. We here briefly sum- 
marize the procedure. All objects were selected 
from the Sloan Digital Sky Survey Data Release 7 
(SDSS DR7) archive according to the following cri- 
teria: (1) redshift in either the 0.35 < z < 0.37 bin 
or the 0.56 < z < 0.58 bin, (2) H/3 equivalent width 
and Gaussian width > 5A in the rest frame. Ob- 
jects with z ~ 0.36 were selected to extend the BH 
mass scaling relations presented in Paper I, II, and 
III to the low-mass range. They meet the addi- 
tional criterion (3) Mbh ^ 1O 8 M as determined 
from the width of the H/3 line and the AL5100 lu- 
minosity measured from the SPSS spectra and 
assuming the calibration given by iMcGill et al 
(|2008l ). 



For two objects in Paper II (0107 and 1015) 
the ACS images revealed dust lanes. Thus, Near 
Infrared Camera and Multi-Object Spectrometer 
(NICMOS) images were additionally obtained to 
correct for extinction. Table Q] summarizes the 
sample properties of all 23 objects. 

3. OBSERVATIONS AND DATA RE- 
DUCTION 

The sample was observed using the NIC 2 cam- 
era and the broad filter F110W (~ J-band) of 
NICMOS on board HST. The 17 objects at z = 
0.36 were observed for a total exposure time of 
2560 seconds per object (11 as part of GO 11208, 
PI Woo; 6 as part of GO 11341, PI Gallagher); the 
6 objects at z = 0.57 were observed for a total ex- 
posure time of 5120 seconds per object (GO 11208; 
PI Treu). Four separate exposures were obtained 
per object, dithering by semi-integer pixel offsets 
to recover resolution lost to under-sampling and 
to improve cosmic-ray and defect removal. 

The individual exposures were first processed 
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with the NICMOS CALNICA pipeline (Version 
4.4.0), and then dither-combined using a custom- 
made pipeline written in IRAF0/STSDAS dither 
package (Version 3.4.2). The pipeline relies on 
the package drizzle and takes care of aligning 
the images, removing sky background, correcting 
for the pedestal effect (the variable quadrant bias 
present in NICMOS images) using "pedsky'H, and 
the NICMOS non-linearity using "rnlincor" , iden- 
tifying and removing cosmic rays and defect pixels 
and finally drizzling all input images together. For 
our final drizzle iteration we chose a drizzle. pixfrac 
parameter of 0.9 and a drizzle. scale parameter of 
0.5, which resulted in a final scale of 0.038 arc- 
sec/pixel. In Figure [TJ postage stamp images of 
all 23 objects are shown. (We refer the reader to 
Paper II for the ACS images of those 2 objects 
in our sample, 0107 and 1015, for which we have 
both ACS and NICMOS images.) 

4. SURFACE PHOTOMETRY 

To decompose nuclear and host-galaxy light 
(spheroid and potentially bar or disk), we used 
GALFIT, a 2D galaxy fitting program that can 
simultaneously fit one or more objects in an im- 
age c hoosing from a library of functio nal forms 
(e.g.. ISersTall968 ; de Vaucouleurslll948l exponen- 
tial, etc.) (jPeng et al.1 120021 ). Decomposition of 
complex images in multiple components is a diffi- 
cult statistical challenge due the degeneracies in- 
volved, and the highly non-linear dependency of 
the likelihood on a large number of parameters. To 
deal with this problem, we develop a methodology 
based on physical assumptions to reduce the num- 
ber of free parameters and extensive trial and error 
exploration of the multidimensional space to avoid 
local minima of the posterior probability. This sec- 
tion describes our fitting procedure in detail. 

4.1. PSF & Noise Image 

For convolution with the point-spread function 
(PSF) of the HST NICMOS optics, and to fit the 



X IRAF (Image Reduction and Analysis Facility) is dis- 
tributed by the National Optical Astronomy Observato- 
ries, which are operated by AURA, Inc., under cooperative 
agreement with the National Science Foundation. 

2 Note that our fields are quite empty and enough blank 
sky is available for an accurate determination of sky and 
pedestal. 



central point source of the AGN, we created PSFs 
using TinyTim (Version 6.3). Compared to other 
cameras onboard HST, TinyTim produces fairly 
good PSFs for NICMOS (and especially NIC2) 
because the PSF is less sensitive to aberrations in 
the infrared ( Krist fc Hook|[2004l ) ; generally, Tiny- 
Tim PSFs are considered an adequate alternative 
when well-matched stellar PSFs are not available 



(|Kim et al.ll2008l ). 

To minimize PSF mismatch due to spatial dis- 
tortion, we simulated PSFs at the location of the 
objects. We created a library of 17 PSFs using 
a wide range of different stellar templates (F6V 
to K7V) and power-law functions (Fx oc A £ with 
e = — 3 to e = 0.5 in increments of 0.5) at the 
four different chip positions of the individual ex- 
posures of the science targets. These four images 
were then dither-combined using the same proce- 
dure as for the science targets. To account for 
breathing, we additionally created PSF models for 
the above range wit h focus values of ±5 /im around 
the nominal focus ([Rhodes et al.ll2007T ). 

We created noise images by dither-combining 
the necessary extension files provided in the image 
block (the output of the CALNICA pipeline; see 
NICMOS data handbook Version 7.0) in the same 
way as the associated science image. 

4.2. Fitting Procedure & Uncertainties 

For each object, we assumed the following 
AGN/host galaxy fitting procedure using GAL- 
FITO We first fitted the central AGN compo- 
nent with a PSF, and thus determined the cen- 
ter of the system, which was subsequently as- 
sumed to be common to all components and fixed. 
We then modeled the sp heroid component with 
a de Vaucouleursl ( 19481 ) profile. We carefully 



checked the images and the residuals for evidence 
of a disk component and added an exponential 
disk if required by the images, residuals, and the 
X 2 statistics. The same approach was used to 
determine the need for an additional bar compo- 
nent, but unlike in Paper II (where seven out of 
17 objects required the fitting of a bar), we did 
not find evidence for a bar in any of the objects 
in our sample. Neighboring objects were fitted si- 
multaneously. Note that the sky was determined 



3 Note that this procedure is in agreement with the one used 
in Paper II and the comparison is therefore straightforward. 
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independently and subtracted out during pipeline 
reduction (see above) which is preferable when 
using GALFIT, as the sky background is not only 
degenerate with the extended wing of the galaxy, 
but it might also be used by GALFIT to com- 
pensate a mis match between intrinsic and fitte d 
galaxy profile (|Peng et alJl2002t iKim et al.ll2008h . 

To ensure that the best resulting fit indeed cor- 
responds to the true global minimum of the \ 2 
over the parameter space, we performed a care- 
ful inter-comparison between different fits with a 
variety of initial parameters and combinations of 
components. Finally, each galaxy was fitted with 
all PSFs in the library to find the best-fitting PSF 
and thus, the best-fitting parameters, and to es- 
timate uncertainties due to PSF mismatch. The 
differences in derived spheroid and PSF magni- 
tudes using different PSFs in our library are small 
(< 0.05 mag) and negligible compared to other 
systematic errors. To understand these systemic 
errors, we simulated artificial images spanning the 
same parameter space as our objects and tested 
how reliably GALFIT can retrieve the different 
components (see Appendix [A]). We use the results 
to estimate our uncertainties and adopt a conser- 
vative total uncertainty on the spheroid luminosity 
of 0.5 mag (i.e., 0.2 dex). The AGN luminosity is 
uncertain to within 0.2 mag. The dominating er- 
ror when constructing the BH mass - spheroid lu- 
minosity relation is the uncertainty of BH masses 
from single-epoch measurements, ~0.4 dex. 

A PSF+spheroid decomposition gives a satis- 
factory fit to the host galaxies of ten out of the 
23 objects (see Table [2]). The remaining 13 ob- 
jects show evidence for a disk component in both 
the residual image and in the \ 2 statistics and 
thus, an additional exponential disk component 
was added. However, in ten of the 13 objects, the 
addition of a disk component results in a vanish- 
ingly small spheroidal component. In these cases, 
we fixed the spheroidal half-light radius to the 
minimum resolvable size of 2.5 pixels (^0.1"), as 
determined by simulations, and consider the mea- 
sured spheroid luminosity an upper limit. For one 
object (1501), even fitting a single spheroid com- 
ponent had the same effect Below, we discuss 



how we combined these upper limits with infor- 
mative priors on the sphcroid-to-total luminosity 
ratio for galaxies to estimate spheroid luminosi- 
ties. 

Note that compared to the ACS images studied 
in Paper II, NICMOS images are less sensitive to 
a potential disk component, dominated by young, 
blue starsQ We cannot exclude to have missed a 
disk component in some of the objects for which 
we only fit a spheroidal component. However, this 
is a conservative approach, i.e. reducing any po- 
tential offset in the Mbh-£sp1i relation. This is 
true in general for our procedure: We only fit a 
two component model consisting of disk+bulge, if 
there is irrefutable evidence for a disk component 
(see e.g. Fig. [IJ. Without such evidence, using 
only one component is conservative in the same 
sense above. 

In Figs. and [31 we show surface-brightness 
profiles derived using the IRAF program "ellipse" , 
for the data as well as each component that was 
fitted using GALFIT. As the fitting was done in 
two dimensions, these figures should be considered 
as illustrations only, showing the relative contribu- 
tions of the different components to the total fit as 
a function of radius. We divide the sample in two 
groups, based on the quality of the fit: In Fig. [2l 
we show all objects that were fitted by a resolved 
spheroid component; in Fig. [31 we show objects 
with an unresolved spheroid component, i.e. those 
for which GALFIT ran into the size limit of the 
spheroid of 2.5 pixels. 

In Appendix [Bl we discuss the effects of the 
choice of different Ser s ic ind ices other than n—4 
(i.e. a Ide Vaucouleursl ( 1948f ) profile). To briefly 
summarize, all results stated in the paper remain 
the same within the errors, when choosing the 
best-fitting Sersic index instead of n=4. 

4.3. Estimating Spheroid Luminosities 
from Upper Limits Using Informative 
Priors 

As described above, for ten objects the addi- 
tion of a disk component resulted in a vanish- 
ingly small spheroidal component. For one object 
(1501), even a single spheroid component had the 



4 Note, however, that fitting this source is particularly com- 
plicated as it is in the process of merging with a neighboring 
galaxy. 



5 At the same time, NICMOS images have the advantage of 
being less affected by dust. 
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same effect @ Thus, for these 11 objects, we fixed 
the spheroidal half-light radius to the minimum 
resolvable size of 2.5 pixels and inferred an upper 
limit to the bulge luminosity. The same is true 
for five objects in Paper II and we thus include 
them in the analysis described here. In brief, we 
combine the upper limit on spheroid magnitude 
with prior knowledge on bulge-to-total luminos- 
ity ratios (B/T) as a function of total host-galaxy 
magnitude. In terms of Bayes' Theorem, we de- 
rive the posterior on B/T by combining our like- 
lihood - in the form of a step function limited to 
the measured upper limit from GALFIT - with a 
prior taken from the literature. 

The prior is determined from quantitative mea- 
surements of the distribution of sphe roid-to-total 
luminosity ratios (jBenson et al.ll2007T ). derived for 
a sample of ~8800 galaxies from SPSS, using 



the 2 D fitting code GALACTICA (jBenson et al 



2002). The galaxy redshifts span a range of 0.02 
< 0.3 with an average of 0.09. The abso- 



< z 



20.2 < i? mag < 



lute R magnitudes of these g alaxies are compara - 
ble to our sample (80% of the Benson et al. (2002) 
galaxies are within our range of 
-22.6)0 

For each object in our sample, we performed 
the follow ing steps. Fi r st, we selected only ob- 



jects from iBenson et al.l (|20071 ). for which the to- 
tal galaxy magnitude is within ±0.5 mag of the 
Seyfert total host-galaxy magnitude (which is typ- 
ically accurate to < 0.1 mag). Then, we created a 
histogram over the B /T values of the Benson et al. 
galaxies within this magnitude range (step size of 
0.1) and cut this histogram at the upper limit of 
B/T we derived for the Seyfert galaxies. Finally, 
we calculated the mean and sigma of the remain- 
ing B/T values and used this value to derive the 
spheroid luminosity for the Seyfert galaxies. In 
Fig. 21 we show the prior, likelihood, and poste- 
rior distribution functions for all 16 objects. The 
upper limit and mean posterior B/T ratio is also 
shown. 



6 As pointed out before, fitting this source is particularly 
complicated as it is in the process of merging with a neigh- 
boring galaxy. 

7 Note that we do not correct for any evolution in luminosity 
here, as the effect is within the errors. 



4.4. Dust Correction 

For two objects (0107 and 1015), the ACS 
F775W images studied in Paper II reveal dust 
lanes in the host galaxy, preventing an accurate 
measurement of the AGN and spheroid luminos- 
ity from the ACS images alone. Thus, NICMOS 
F110W images were obtained, which we use here 
for dust correction. Briefly, the color excess is 
measured from the two colors and used to cor- 
rect for dust extinction assuming an extinction 
law. The proce dure we adapt is simil ar to the 



one described in lKoopmans et al ( 2003 ) 



First, we deconvolved the ACS image using the 
ACS PSF from Paper II and a Lucy-Richardson 
algorithm (IRAF program "lucy"). Then, the 
F775W image was rotated to match the orienta- 
tion of the NICMOS image and drizzled to the 
same pixel scale. We assume that the centroids 
are unaffected by the dust lane - an assumption 
supported by the distribution of the dust seen in 
the images - and we thus centered both images 
on their peaks. In the next step, the F775W im- 
age was convolved by the NICMOS PSF (IRAF 
program "imconvolve" ) to match the resolution of 
the F110W image. Then, a color map was created 
from the ratio of these matched images and the 
intrinsic color was assumed to correspond to the 
minimum of the color and to be spatially uniform. 
Finally, the color excess was converted into ex- 
tinction assuming Ay = 3.1Eb-v> A-pnow/Ay — 
0.628 (corresponding to rest-frame F814W), and 
Afyysw/Av — 1.049 (corresponding to rest-frame 
F555W). The extinction-corrected NICMOS im- 
age was used for fitting with GALFIT. 

5. DERIVED QUANTITIES 

5.1. Rest-Frame V-Band Luminosities 

We applied correction for Galactic extinc- 
tion, assumin g Ay = 3.1E(p_y) and Apnow = 
0.902E (B _ V) (ISchleeel et al.lll998l). The values for 



Eb-v were taken from lSchlegel et al.l (|1998I ). The 



F110W AB magnitudes were transformed to rest- 
frame optical bands by performing synthetic pho- 
tometry on an early-type galaxy template spec- 
trum, a procedure traditionally referred to as k- 
correction. The template spectrum initially has 
arbitrary units, and these units were adjusted so 
that the synthetic observed frame F110W magni- 
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tudes match the magnitudes from our photometry. 
We then evaluated the V-band magnitudes at the 
rest-frame of the template; luminosities were de- 
termined by correcting for the distance modulus 
given our adopted cosmology. The errors on ex- 
tinction and rest-frame transformation are a few 
hundredths of a magnitude. We note that the 
F110W band roughly corresponds to the R and I 
bands for our two intermediate-redshift samples; 
considering the small scatter in the red colors of 
bulges (that is, the V-R and V-I colors) we are 
able to determine robust estimates of the V-band 
magnitude. We estimate an uncertainty of <0.05 
mag (using the scatter in 20 single stellar popula- 
tion templates with ages ranging from 2 Gyr to 8 
Gyr). 

5.2. Luminosity Evolution 

To allow a direct comparison of the observed re- 
lation in the more distant universe and local sam- 
ples, we evolved the spheroid luminosity according 
to the evolution measured from the evolution of 



the fundamental plane bv lTreu et al.l (|20Qll ): 



logiy :0 = logLy - (0.62 ±0.08 ±0.04) x z (1) 

This corrects pure passive luminosity evolution, 
i.e. the decrease in spheroid luminosity due to an 
aging stellar population. We used the same cor- 
rection for our intermediate-z Seyfert sample and 
for the local RM AGNs which will be discussed 
in § 16.11 However, equation [1] is only valid below 
z < 1 and an equivalent measurement is not avail- 
able at higher redshift. Therefore, for the high-z 
sample we adopt a conservative correction based 
on the predicted evolution for a maximally old 
stellar population (see § 16. 2|) . 

5.3. Black Hole Mass 

As in Paper I, II, and III of this series, black 
hole masses were estimated using the empiri- 
cally calibrated photo-ionization method (e.g., 



Wandel et al.ll999HVestergaardl2002|IWoo fc Urry 
2002f IVestergaard fc Peterson! 120061 : iMcGill et aTl 



for a detailed discussion see Paper II). 
Briefly, the method (sometimes called the "virial" 
method) assumes that the kinematics of the 
gaseous region in the immediate vicinity of the 
BH, the broad line region (BLR), traces the grav- 
itational field of the BH. The width of the broad 



emission lines (e.g. H/3) gives the velocity scale, 
while the BLR size is given by the continuum 
luminosity through application of an empirical re- 
lation found from reverberation ma pping (RM 



e.g., Wandel et al. 19991 : Kaspi et "aTl 



(RM) 
2005; 



Bentz et al.ll2006l) . Combining size and velocity 



gives the BH mass, assuming a dimensionless co- 
efficient of order unity to describe the geometry 
and kinematics of the BLR (sometimes known as 
the "virial" coefficient). This coefficient can be 
obtained by matching the Mbh-c relation of local 
active ga lactic nuclei (AGNs) to that of quiescen t 
galaxies (|Onken et all 120041 : iGreene fc Holl2006h . 
Alternatively, the coefficient can be postulated 
under specific assumptions on the geometry and 
kinematics of the BLR. We note that the exact 
value of the virial factor does not affect our re- 
sults since the relative offset between local and 
higher redshift AGNs should be independent of 
the virial factor. 

We use the following formula which includes 
calibrations of the BLR size-luminosity rela- 
ti on (after subt raction of host galaxy light; 



from 



Bentz et al.1 2006) and a virial coefficient taken 



Onken et alj (2004): 



log -Mbh = 8.58±21og 



3000 km s" 



-0.518 log 



XL 



5100 



10 44 ergs" 



(2) 

with cth/3 the second moment of the broad H/9 
emission line and AL5100 the observed nuclear lu- 
minosity at 5100A. 

To obtain the 5 100 A continuum luminosity 
of the AGN, we extrapolated the extinction- 
corrected PSF AB magnitude in F110W to rest- 
frame 5100A, assuming the power law f^ocv . 
The power-law index -0.5±0.15 is the average 
derived for the majority (24/40) of our sample 
from SDSS photometry in the restframe wave- 
length range 5000A-6600A (as covered by SDSS 
for z=0.36) after subtraction of host-galaxy light 
contribution (Szathmary et al., in preparation). 
Within the error s, this slope is in agreeme nt with 
other studies (see lVanden Berk et al. 2001 and ref- 
erences therein). The uncertainties on the slope 
lead to an uncertainty of <10% in L5100 and <5% 
in Mbhj i-e. negligible compared to the uncer- 
tainty of 0.4 dex when deriving BH masses from 
single-epoch spectra. 

The line width anp was derived from spec- 
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tra obtained with the Keck telescope, using the 
longslit spectrograph LRIS to measure the stellar 
velocity dispersion (see Paper I and Woo et al. 
2009, in preparation, for details). We assume a 
nominal uncertainty of the BH masses measured 
from single-epoch spectra of 0.4 dex. 

Note that we do not corre ct for possible effects 



of ra diation pressu re (e.g. , Marconi et al. 2008, 



2009, see, however, |NetzerJ [2009j) . First, the role 



of radiation pressure on the measurement of BH 
masses is still discussed controversially and sec- 
ond, neglecting its effects is a conservative ap- 
proach: If radiation pressure does affect the mo- 
tion of the BLR clouds, not taking it into account 
would lead to an underestimation of the BH mass. 
Thus, including radiation pressure, the observed 
offset would further increase. 

All results are summarized in Table [2] In addi- 
tion to the 23 objects in the sample studied here, 
we give the results for the sample in Paper II, 
which changed slightly (<0.15 mag) due to a small 
error in the extinction correction in Paper II (con- 
sequently, also the derived BH masses changed 
slightly). Also, we determined spheroid luminosi- 
ties for those objects in Paper II that only had 
upper limits by applying the informative prior. 
Moreover, Paper II used the B-band luminosity of 
the spheroid component (for comparison with in- 
active galaxy samples in the local Universe). We 
here give the V-band luminosity of the spheroid 
component. 

6. COMPARISON SAMPLES 
6.1. Local Comparison Sample 

Interpreting the AfBH-L sp h relation for the dis- 
tant Seyfert samples studied here and any possi- 
ble evolution with redshift requires a robust local 
baseline - ideally of Seyfert galaxies with com- 
parable BH masses and spheroid luminosities to 
avoid selection biases as much as possible. The 
most appropriate local comparison sample for our 
study is the reverberation-mapped sample of 35 
AGN hosts. This sample has the great advantage 
that the BH mass is derived directly via RM and 
does not depend on the BLR size-luminosity re- 
lation and its uncertainties. A detailed analysis 
of HST images of the RM sample to derive AGN 
an d spheroid luminosit ies was recently completed 
bv iBentz et al. ( 2009bl ). However, a comparison 



with our study is not straightforward due to their 
very different approach which aims to get the best 
estimate of AGN-to-host-galaxy luminosity ratio. 
In particular, while we use the simplest decom- 
positi on possible, i.e. PSF, s pheroid (modeled as 
simple de Vaucouleursl (Il948l) p r ofile) p lus possible 
exponential disk, IBentz et al.l (|2009bl ) allow the 
Sersic index to vary, sometimes include more than 
one PSF for the same galaxy (to compensate for 
PSF mismatch), and up to three different spheroid 
components. 

We thus decided to perform an independent 
analysis of the HST archival images, for a homo- 
geneous treatment of all the data, using the same 
approach as for our distant Seyfert galaxies. Our 
robust subsample of the RM AGNs consists of 19 
objects (10 PG quasars and 9 Seyfert-1 galaxies). 
In Appendix [Cj we summarize the details of the 
analysis, and show, 



that our results are in over- 
r 



all agreement with those of IBentz et al.l (|2009bl ). 
However, the spheroid luminosities we der ive are 
often brighter than in IBentz et al. I (|2009bh . espe- 



cially in those cases where we fit a sph eroid compo 
nent on ly and not spheroid+disk as in lBentz et al 



( 2009bl ). since we did not find evidence for a disk 
component. 

As for the more distant sample, we correct 
for passive luminosity evolution to zero redshift 
(§ 15. 2p . This is important as the most massive 
black holes in the RM sample are systematically 
found at higher redshift (up to z=0.29, look-back 
time ^3.3 Gyrs), which changes the b est fit M rh- 



L sn h wi th respect to that presented bv lBentz et al 



2009J3), who did not take luminosity evolution 
into account. 

To comp a re our local relation to that of 
Bentz et al.l (|2009al) . we fitted the data usin g 
the BCES algorithm (jAkritas &: Bershadvl fl996h . 
which takes into account the effects of errors on 
both coor dinates using bivari ate correlated errors. 
Following IBentz et al.l ( 2009a ). we adopt the boot- 
strap of the BCES bisector value with N = 1000 
iterations. We give the different fits in Table [4] in 
the form of 



log 



Mi 



BH 



10 8 M,r 



K + a log 



L. 



sph,V 



10 10 L 



(3) 







However, in line with the Bayesian approach fol- 
lowed in this paper, instead of using the BCES fit- 
ting routine to determine our standard local base- 
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line, we apply our own fitting routine. Following 
standard procedures, gaussian errors on both vari- 
ables are taken into account. The intrinsic scatter 
is a free parameter and is modeled as a gaussian 
distribution. Uniform priors are assumed on each 
free parameter. The inferred slope and intercept 
after marginalizing over the intrinsic scatter are 
given in Table [U As can be seen, the resulting 
slope can range from a = 0.67 to a = 0.81, de- 
pending on the evolutionary correction and on the 
fitting technique. 

For comparison, we transformed the B-band 
magnitu des of the local in a ctive c omparison sam- 
ple from iMarconi k. Hunt! ( 20031 ) (group 1 only) 
to V-band (as suming an elliptical te mplate and B- 
V=0.96 mag; iFukugita et al. 19951) . Using again 
our linear fitting routine including gaussian er- 
rors and intrinsic scatter gives a steeper slope of 
a=l.ll±0.13 (K=0.07±0.08; scatter=0.38±0.07). 
For a discussion of the difference in slope be- 
twe en AGN sample an d inactive galaxy sample, 



sec 



Bentz et al.l (2009a). 



6.2. High-z Comparison Sample 

To study the evolution of the AfBH-^sph rela- 
tion, we selected a h i gh-red shift comparison sam- 
ple from lPeng et a l. (2006bj), consisting of a total 
of 31 gravitationally lensed quasars and 20 non- 
lensed quasars at redshifts of 0.66 < z < 4.5. It in- 
cludes 15 non- lensed (radio-loud an d radio-quiet 



quasars taken from the li terature (jKukula et al 
20011 iRidgwav et"aDl200ll ). We exclude four ob- 
.jects which were al so excluded from the analysis in 
Peng et al.l (|2006q ). one object, for which the BH 
mass is only a lower limil|ff (PSS 2322+1944), as 
well as two extreme outliers in Mbh with high un- 
certainties (B2045+265 and HE 2149-2745), leav- 
ing us with a sample of 44 QSOs (17 non lensed 
and 27 lensed objects). BH masses were estimated 
from single epoch spectra using the broad lines 
C IV, Mg II, or H/3. The luminosity was de- 
rived from two-dimensional su rface brightness fit- 
ting (GALFIT). Note that the lPeng et all (|2006bh 
measurements comprise the total host-galaxy lu- 
minosity instead of spheroid luminosity alone, as 
only one component was fitted. However, there 
is no evidence for any of the objects to have two 
components, indicating that the host galaxies are 



ellipticals. (Note that even if there was a disk 
component present, fitting only one component is 
a conservative approach in the sense that the offset 
from the local relation is the smallest.) 

To allow for a homogeneous treatment of the 
data, we corrected the BH mass estimation based 
on the Mg II line for normaliza t ion di fferences us- 
ing the recipe bv lMcGill et al l d2008l) . Note that 
for H/3 and C IV, IPeng et al.l (|2006bh used a com- 
parable normalization factor and the difference is 
negligible. For objects for which both C IV and 
Mg II measurements are available, we use the lat- 
ter line, as determining Mbh based on the former 
may have larger uncertainties: The C IV line is of- 
ten found to be blueshifted with a strong blue ex- 
cess asy mmetry indicating an outflow component 
(see e.g. Baskin fc Laorlliooih . 



To correct for luminosity evolution, but lacking 
direct determinations of passive evolution out to 
these redshifts, we apply a conservative evolution 
correction based on maximally old stellar popula- 
tions. Specifically, we assume that the single burst 
occurred at z = 5. We use a Salpeter initial mass 
function with solar metallicit y and stellar popu- 
lation s synthesis models from iBruzual &: Chariot 



( 2003 ) to compute the evolutionary correction^ 



tion which might be present (e.g. Martel et al. 




1999 


; Canalizo & Stocktonl 12001; Evans et al. 


2001 


; Scoville et all 120031; iKauffmann et al. 20031 




Sanchez et al.l2004 


Tadhunter et al.|2005t Canalizo et al. 


2006; Barthef 2006 


; Jahnke et al. 2007) - for ex- 



ample due to triggering of SF from a merger event 
that also triggered the AGN activity - we would 
infer faster passive evolution and therefore more 
pronounced evolution in Mbh at fixed present-day 
luminosity. Thus, our approach is conservative. 

7. RESULTS 

7.1. Host-Galaxy Morphology And Merger 
Rates 

From the final reduced images (see Fig. [J), we 
derive the overall host-galaxy morphology. At 
least five of the 23 objects show a clear large spi- 



3 Assuming Eddington-limited accretion 



9 We used an Sbc template for obtaining V-band rest-frame 
magnitudes (instead of an elliptical template used for the 
intermediate- and low-z sample), as it is closer to the colors 
predicted by our model. 



9 



ral disk (0804, 1043, 1046, 1352, 2340) and one 
object has an extended disk-like structure (1007). 
Including the sample from Paper II, a significant 
fraction of the host galaxies (> 15/40) have mor- 
phologies of Sa or later. 

In the NICMOS images, seven objects show ev- 
idence for tidal interactions and merging such as 
tidal tails and other morphological disturbances 
(0934, 1021, 2158, 0155, 0342) or nearby com- 
panions connected by tidal features (1501, 1526). 
Half (3/6) of the objects at z=0.57 show signs of 
mergers/interactions (0155, 0342, 1526). In some 
cases, it is difficult to clearly distinguish between 
the presence of a spiral disk and tidal tails; we 
cannot exclude the presence of a tidally disrupted 
disk (e.g. 1352, 0342). Combined with the objects 
in Paper II, 13 of 40 objects show some sign of 
tidal disturbance. When considering only objects 
at z=0.36, 10/34 objects are in apparently dis- 
turbed systems (0.29±0.1). This agrees with the 
fraction of disturbed systems found for a control 
sample of GOODS galaxies: Selecting all galaxies 
within GOODS with comparable stellar luminos- 
ity and a redshift range of z = 0.36 ± 0.1 and 
performing the same visual classification lead to 
12/42 disturbed systems (0.28 ± 0.08; Paper II). 
This fraction is somewha t larger than obse rved 
in the lo cal universe (e.g ., IPatton et al. 20021 . see, 
however, 



Tal et al. 200' 



7.2. BH Mass - Spheroid Luminosity Re- 
lation 

The resulting MBH-^sph relation is shown in 
Fig. O (upper left panel). Objects with signs of 
tidal interaction or merger are marked by open 
black circles; they are not significant outliers. In 
the upper middle panel of Fig. [5j the high-z sam- 
ple is included. Fig. [5] (upper right panel) shows 
the distribution of the residuals in logAfsH with 
respect to the fiducial local relation. If we treat 
the intrinsic scatter of the relation as a free pa- 
rameter and marginalize over it, the offset we de- 
rive with respect to the local relation (solid black 
lines in Fig. El upper left panel) is A log Mbh = 
0.14 ± 0.07 ± 0.20 (statistical and systematic er- 
ror; w.r.t. L sp h,v including the full sample at both 
z = 0.36 and z = 0.57). For comparison, in Pa- 
per II we found A log Mbh = 0.51 ± 0.14 ± 0.19 
(i.e. when considering only the blue data points 
in Fig. upper left panel). Although the num- 



bers are consistent within the errors, we note that 
they cannot be compared directly due to the rad- 
ically different selection function in Mbh, for the 
two samples (see § 17. 3ft . Expressed as offset in 
spheroid luminosity, AlogL S ph = -0.19 ± 0.08 ± 
0.21. 

Studying the evolution with redshift of the off- 
set in BH mass, with respect to the fiducial lo- 
cal relation, we fit a linear relation of the form 
AlogA/sH = 7log(l + z) and include the intrin- 
sic scatter in log A/bh as a free parameter. We 
assume negligible errors on the redshifts and a 
standard error of 0.4 dex for log Mbh- We find 
7 = 1.3 ± 0.9. Note that this fit does not take 
into account systematic errors nor selection ef- 
fects. Ad ding the higher red shift comparison sam- 



ple from iPeng et al. I (|2006blh we extend our red- 



shift baseline, decreasing the error on the slope, 
resulting in 7 = 1.2 ± 0.2 (Fig. H left panel). We 
also plot the evolution as a function of look-back 
time (Fig.[6l right panel). However, we stress that 
the figures and the fits discussed in this section 
are for illustration only as they ignore selection ef- 
fects. The correct quantitative results taking into 
account selection effects are given in § 17.31 

Note that without correction for passive lumi- 
nosity evolution, there is little to no offset for any 
of the distant objects. A s i milar r esult has already 
been found by Peng et al. ( 2006bl ) , who show that 
host galaxies harboring BHs of the same mass were 
as luminous at a redshift out to z=4 .5 as they are 
today , up to ~12 Gyrs later (see also lDecarli et al 
2009h . 



7.3. Selection Effects 

As discussed in Paper II, by selecting targets 
based on their nuclear properties and in particu- 
lar on the presence of a broad line A GN, we may 
be bi asing our inferred offset fs ee alsolLauer et alJ 
l2007t ). an effect analogous to the lMalmquistl (|l924|) 
bias. The magnitude and sign of the bias depends 
on the errors, on the selection function, on the 
spheroid luminosity function and on the intrinsic 
scatter of the correlation with host-galaxy lumi- 
nosity. Here, we exploit the larger sample to cor- 
rect the offset and infer for the first time the intrin- 
sic scatter of the relation. The slope is assumed to 
be fixed to the local value, given that the dynamic 
range of the data is not sufficient for an indepen- 
dent determination of its evolution. 
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Briefly, we use a Monte Carlo approach to sim- 
ulate the observations including selection effects 
and compute the likelihood and posterior distri- 
bution function as a function of the two free pa- 
rameters: slope (3 of the relation A log Mbh= 
(3 log(l + z) at fixed zero redshift spheroid luminos- 
ity, and intrinsic scatter cr int of the MeH-isph re- 
lation which is assumed to be non evolving. First, 
we populate the local MBH-L sp h correlation ac- 
cordi ng to the sph e roid luminosity function taken 
from Driver et al. (|2007l) (their Table 1, Sample 
"Ellipticals + bulges"). Second, for each value 
of the free parameters (3 and tJi n t, and for each 
object in the distant sample, we generate a sim- 
ulated observed sample, assuming gaussian errors 
on both axes, with amplitude equal to the observa- 
tional errors. Third, we model the selection effect 
by hard thresholds in log Mbh, as appropriate for 
each sample: For the initial samples of z = 0.36 
Seyferts introduced in Papers I and II, as well as 
for the sample at z = 0.57, we adopt the inter- 
val [7.5,9] (lower and upper value of log Mbh); for 
the z = 0.36 sample introduced here - which was 
selected to have small black hole masses and was 
thus restricted to log Mbh ^8.2 - we use [7.5,8.2]; 
for the high-z sample we assume [7.5,10]. It is im- 
portant to notice that both the upper and lower 
limits are relevant for the analysis, as they bias 
the results in opposite directions. For each ob- 
ject in the distant sample, we select simulated ob- 
jects with consistent spheroid luminosity within 
the error, generate a one dimensional simulated 
distribution in Mbh, and compute the likelihood. 
To compute the posterior, we adopt a uniform 
prior on (3 and two choices of prior for cr^: i) 
uniform, appropriate when the parameter is un- 
known but believed to be of orde r unity; and ii) 
aint = 0.38±0.09, as determined bv lGultekin et al 



( 20091 ) for a local sample of inactive galaxies. The 
latter is the most informative prior, although it 
comes with the price of assuming that the scatter 
of the MBH-^sph relation is the same for active and 
inactive galaxies. The results are shown in Fig. [7] 
If the high-z sample is included, (3 — 1.4 ± 0.2 is 
well determined regardless of the assumed prior 
on (Ti n t. For a uniform prior on Cinti the inferred 
scatter is 0.3±0.1 dex (<0.6 dex at 95% CL). If 
the high-z sample is not included, the baseline in 
redshift is not sufficient to determine [3 and <7i n t 
simultaneously for a uniform prior. Adopting the 



prior by iGultekin et al. (|2009h to break the de- 
generacy, we find /3=2.8±1.2. This is consistent 
with the trend observed for the complete sample 
although the bounds are weaker due to the smaller 
baseline in redshift. 

When excluding those 16 objects in our sam- 
ple for which we estimated spheroid luminosities 
from upper limits using informative priors (§ 14. 3| . 
we obtain /3=1.3±0.3 (including the high-z sam- 
ple) and /3=1.8±1.4 (without the high-z sample 
using the same prior as above), i.e. consistent 
within the errors. If we exclude all objects for 
which Mbh was estimated based on the C IV line 
(which may be more uncertain; § 16. 2|) . the evolu- 
tion is less well constrained, since half of the high- 
z objects are excluded. Using again a prior on 
(Tint = 0.38±0.09 as above, it results in /3=1.1±0.3. 
The slope also gets shallower w hen using the lo- 
cal in active galaxy sample from Marconi fc Hunt 
(2003) ( group 1 only, transformed to V-band mag- 
nitude, see §[63]): /3=0.9±0.2. 



7.4. 



BH Mass - Host- Galaxy Luminosity 
Relation 



We calculate the total host-galaxy luminosity 
for both our intermediate-redshift Seyfert sam- 
ple and the local RM AGNs and show the Mbh- 
ihost relation in Fig. [5] (lower left panel). Note 
that, for consistency and lack of additional infor- 
mation, we assume the same k-correction template 
and passive luminosity evolution for the total host 
galaxy as for the spheroid luminosity (see § 15. 2| . 
Conservatively, we also assume the same error on 
the total luminosity as on the spheroid luminosity 
of 0.5 mag, although, generally, the error on the 
total luminosity is smaller. 

Compared to the Mbh-^spIi relation, the A/bh- 
ihost relation is apparently non-evolving: If we 
again treat the intrinsic scatter of the relation as 
a free parameter and marginalize over it, the offset 
we derive with respect to the local relation (solid 
black lines in Fig.[H lower left panel) is A log A/bh 
= -0.03 ± 0.09 ± 0.04 (w.r.t. L hos t,v; including the 
full sample at both z = 0.36 and z — 0.57). Ex- 
pressed as offset in spheroid luminosity, A log L sp h 
= 0.04 ± 0.09 ± 0.04. 

The best fit to the local RM AGNs (black solid 
line in Fig. [5j lower left panel) gives a marginally 
steeper slope than for the MBH-^sph relation (a — 
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0.96±0.18 vs a = 0.70±0.10; Table gj). Overplot- 
ting the high-z comparison sample (Fig. [5] lower 
middle panel), their luminosity remains the same 
as in t he upper middle pan el: The objects were fit- 
ted by P eng et al. (200Gb) by only one component 
(without any evidence of a second component) and 
thus L sp h=£host- Apparently, the high-z compari- 
son sample does not follow the same A/BH-^host re- 
lation, instead the offset remains. The distribution 
of the residuals in logA/eH of the distant AGNs 
with respect to this fiducial local relation is shown 
in Fig. [5] (lower right panel). 

8. DISCUSSION 

8.1. The role of mergers 

Theoretical studies generally invoke mergers 
to explain the observed scaling relations between 
BH mass and host-galaxy spheroid properties 
- a promising way to grow both spheroid and 
BH. In a simple scenario, spheroids grow by (i) 
the merging of the progenitor bulges (assuming 
that both progenitors have a spheroidal com- 
ponent), (ii) merger-triggered starbursts in the 
cold galactic disk, and (iii) by tran sforming stel- 
lar disks into stellar sphe r oids (e.g.. |Barneslll992t 



Mihos fe Hernquistl Il994t ICox et all 12004 ). thus 
increasing the spheroid luminosity and stellar ve- 
locity dispersion. The fueling of the BH, on the 
other hand, is triggered by the merger event as 
the gas loses angular momentum, spirals inward 
and eventually gets accreted onto the BH, giving 
rise to the bright AGN o r 'quasar' period in the 
evolution of galaxies (e . g^, iKauffmann fc Haehnelt 



20001 : iDi Matteo et al.ll2005h . Eventually, if BHs 
are present in the center of both progenitor galax- 
ies, they may coalesce. In such a simple scenario, 
an evolution in the BH mass - spheroid luminosity 
relation is not necessarily expected: Both spheroid 
and BH grow from the same gas reservoir, and 
bulge stars added to the final spheroid followed 
the BH mass - spheroid luminosity relation prior 
to merging, so the relation will be preserved when 
the BHs coalesce. However, while mergers pro- 
vide a way to grow both spheroids and BHs, they 
may do so on very different timcscales. Moreover, 
the merger history of galaxies varies, depending 
e.g. on formation time and environment. Different 
types of merger, for example with a different rela - 
tive role of dissipation (e.g., Hopkins et al. 20091 ) 



have different effects on the growth of spheroid 
and BH: For a gas-rich major merger between an 
elliptical galaxy and a spiral galaxy - the latter 
without a (massive) BH - the bulge grows more 
efficiently t han the BH b y the disruption of the 
stellar disk (|Crotonll2006h . 

In general, our images of the intermediate-z 
Seyfert galaxies support the merger scenario (see 
Fig. [1] and § 17. ip . However, objects with evidence 
for merger/interaction do not form any particular 
outliers in the BH mass - spheroid luminosity re- 
lation (Fig. [5j. This may not be too surprising: 
For those objects for which we still see two sepa- 
rate galaxies in the process of merging, we fitted 
both separately and the bulge luminosity of the 
AGN host has not yet increased from the process 
of merging. Other objects with signs of interac- 
tion may be in a later evolutionary stage where the 
bulge luminosity has already increased and thus, 
the object falls closer to the local relation. Finally, 
mergers between similar objects would only move 
the system parallel to the local relation. In gen- 
eral, the effect of mergers on the measured bulge 
luminosity of an object depends on the type of the 
merger, the evolutionary stage of the merger, and 
the timescales involved to grow spheroid and BH. 
Such a detailed comparison of merger type and age 
with theoretical predictions is beyond the scope of 
this paper, given the small sample of merging ob- 
jects and the limited information at hand. 

Note that the fraction of apparently disturbed 
systems we find is not higher than that of a 
comparison sample of inactive galaxies at the 
same redshift ( fl7.ip . Thus, from our images 
alone, we cannot infer a causal link between 
a merger/interaction event and the AGN ac- 
tivity we observe. Instead, "normal" galax- 
ies may have the same merger history, and on- 
going interactions arc not necessarily predic- 
tive of AGN activity. The role of mergers for 
the fueling of AGNs is debat ed in the liter- 
ature (e.g. ISanders et al.1 1988 



1984; Hutchings et al 



Hcckman et al 
Disney et al. 1995; 



Bahcall et al.lll997l:]McLure et al.ll999tlCanalizo fc Stockton 
200l [ lDunlop et al.M2003tlFloyd et al.l2004|CanaIizo et all 



20071: lUrrutia et all 120081: iBennert et all |2008; 
Veilleux et all 120091 : iTal et al.ll2009h . While there 
is little doubt that mergers are helpful, they are 
certainly not a sufficient condition, considering the 
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numerous inactive interacting galaxies0 Also, 
mergers may be necessary for the high-luminosity 
QSOs only while for Seyfert galaxies, secular evo- 
lution through processes such as bar instabilities 
may be the dominant effect in the evolution of 
these galaxies. We will come back to this issue in 

31 



8.2. BH Mass - Spheroid Luminosity Re- 
lation 

Combining results of low-z, intermediate-z and 
high-z AGNs, treated in a self-consistent man- 
ner, we can estimate the intrinsic scatter of the 
MBH-^sph scaling relation and correct evolution- 
ary trends for selection effects. We discuss scatter 
and evolution in the next two subsections. 

8.2.1. Scatter of MeH-isph 

The intrinsic scatter we find (0.3±0.1 dex; <0.6 
dex at 95% CL) is non-negligible. However, we as- 
sume the intrinsic scatter of the M^n-Lsph relation 
to be non evolving. While it would be desirable 
to directly study the evolution of the scatter with 
redshift, this requires a larger sample than the 
one we have at hand. Actually, we might ex- 
pect a larger intrinsic scatter at higher redshifts, 
given the different ways and timescales involved 
when growing spheroids and BHs through merg- 
ers. Indeed, for the local Universe, the observed 
tightness in the relations has been a challenge for 
theoretical studies. It has been ex plained by self- 



regula ted models of BH growth (jHopkins et al 
l2009bl) in which the energetic feedback of the AGN 
eventually halts accretion, preventing the BH 





Ciotti & Ostrikerl Il997l 120011 ISilk & Rces 


1998: 


Murrav et all 12005: 


Di Matteo et al. 


2005; 


Sazonov et al.l2005; Hopkins et al.ll2005; Suringel et al. 


2005; Di Matteo et al. 20081). 



Also, a significant fraction of the host galaxies 
of both our local RM AGN sample (~9/19) and 
our intcrmcdiate-z sample (> 15/40) are promi- 
nent late-type spirals of type Sa or later which 
have been found to have a larger intrinsic scat 



ter than elliptical galaxies (e.g., iGiiltekin et al 



10 However, this might also be due to the timescales involved, 
with the signs of interaction outliving the A GN activity 
(see also the case of present-day Type II AGNs; Choi et al. 
e.g. 120091) . 



20091 for the Mbh-c relation: 0.44 dex when in- 
cluding all galaxies vs. 0.31 for elliptical galax- 
ies only). As already discussed in paper II, the 
intermediate-z late-type spirals may eventually 
fall on the local relation later, through merging , 



in line with "downsi z ing" ( e .g. Cowie et al.l 



Brinchmann fc Ellisl l200d: Kod ama et alJ 
Bell et all 120051 iNoeske et al.l l2007h : Less 



1996; 



2004; 



mas- 
sive, blue galaxies merge at later times and arrive 
at the local relation by becoming larger, bulge- 
dominated red galaxies. Also, at least some spi- 
ral galaxies may not have classical bulges, but 
pseudobulges which are characterized by surface- 
brightness profiles closer to exponential profiles, 
ongoing star formation or starbursts, and nuclear 
bars or spirals. It is generally believed that they 
have evolved secularly through dissipative pro- 
cesses rather than being formed by merge rs (see 
e.g. review m iKormendv fc Kennicuttlliool . BHs 
have been found to reside in galaxies without clas- 
sical bulges w hich may not follow the same scaling 
relations (e.g., Greene et al. 2008f h 



8.2.2. Evolution o/MBH-^sph with redshift 

To generalize our results and to facilitate com- 
parison with theoretical and observational works, 
it is useful to estimate the evolution of the Mbh- 
spheroid stellar mass relation. We can convert the 
observed evolution of A/bh - spheroid luminosity 
into that between Mbh and spheroid mass, if we 
assume that - after correction for luminosity evo- 
lution - the mass-to-light ratio does not change 
from sample to sampler H Under this assumption, 
an offset of AMbh at fixed L sp h equals that at 
fixed M Btar and thus, M BH /M sph oc(l + z) 1A±0 - 2 . 

We are now in a position to make a broad 
range of comparisons. In the literature, the BH 
mass evolution is d iscussed quite controversially. 
Shields et all (l2003h study the Mbh-c relation out 
to z — 3.3, estimating A/bh from H/3 and a from 
[O III] and find that the QSOs and their host 
galaxies follow the local relation. (Note, however, 
that using [O III] as a surrogate for a can be prob- 
lematic as [O III] is known to often h ave an outflow 
component; for a discussion se e e.g. I Greene fc Ho 
(|2005l ); lKomossa fc Xul (|2007[h l A similar conclu- 



11 Unfortunately, spatially resolved color information for a 
more sophisticated estimation of the stellar mass of the 
bulge is currently not available. 
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sion has been reached by IShen et all (|2008h who 
study over 900 broad- line AGNs o ut to z ~ 0.4 
from SDSS. lAdelberger fc Steidell (j2005h use the 
correlation length of 79 quasar hosts at z ~ 2 — 3 
to estimate the virial mass of the halo and the C IV 
line width and UV flux at 1350A to estimate Mbh- 
When comparing th e resulting Mrh -Mhaio rela- 
tion to the local one ( Ferraresa 120021 ). thev find no 
evidence for evolution. In particular, they can rule 
out evolution of the form M sp h/M sp h oc(l + z) 2 - 5 
with z = 2.5 at 90% CL, given their error bars. 

Other observational studies find the same 
trend in evolution as we do, i.e. that BHs are 
too massive for a given bulge m ass or veloc- 
ity d is persion at h i gher r e dshifts d Walter et al. 
2004 (shields et all bood iMcLure et all l2006t 



Peng et afT2006bHSalviander et al.ll2007tlWeiss et al 



20071: iRiechers et al.|l2008l |2009HGu et al. Il2009h .~ 

McLure et al.l (|2006l ). for example, study radio- 
loud AGN (0 < z < 2) ari d find M BH /M SDh cx(l + 
z )2.o7±o.76 |Peng et alj (|2006bl ). whose data, 
treated in a consistent manner to match our data 
set, are included in this study, rule out pure lumi- 
nosity evolution and find that the ratio between 
A/bh and M sp h was ~four times larger at z ~ 2 — 3 
than today. For 89 broad-line AGN s between 1 < 



z < 2.2 in the zCOSMOS survey, iMerloni et al 
(|2009t ) find M BH /M sph cx(l + g) - 7 ** - 12 . How . 



ever, this fit refers to the total host galaxy in- 
stead of the spheroid component alone. At least 
some galaxies will have a non-negligible disk frac- 
tion, which, when taken into account, would re - 
sult in a larger offset (see also lJahnke et al1 l2009). 
For a sample of ~ 100 quasars sele c ted to re- 
side in elliptical hosts, Decarli et al. ( 20091 ) es- 



timate that M-qh /Msph was ^8 times larger at 
z ~ 3 than today, i.e. Mbh / M sp h(x(l + z) 1 - 5 . 
The evolution we find is also consistent with our 
previous results, within their much larger errors: 
M BH /M sph cx(l + z) 15±1 from the M BH -L aph re- 
lation in Paper II, whose data are included here, 
and A log M BH oc(l + z ) 31±1 - 5 from the M BH -(j 
relation in Paper III. 

From a theoretical perspective, the discussion 
on the evolution of the BH mass scaling relation s 
is not any less controversial. IShankar et all (|2009l ). 
for example, use the local velocity dispersion func- 
tion (VDF) of spheroids, together with their in- 
ferred age distributions, to predict the VDF at 
higher redshifts. Using the M B h-c relation with 



a normalization allowed to evolve with redshift 
(cx (1 + z) s ), they infer the BH mass density and 
compare it to the accumulated BH mass density 
derived from the time integral of the AGN LF. 
They find a mild redshift evolution (<5 < 0.35), 
excluding 5 > 1.3 at more than 99% CL (with 
the possibility of a stronger evolution for the more 
massive BHs). Another study using fully cosmo- 
logical hydrodynamic simulations of ACDM fol- 
lowing the growth of galaxies and supermassive 
BHs, as well as their associated feedback pro- 
cesses, finds only limited evolution in M BB wit h 



a steepening at z= 2-4 (|Di Matteo et al.l l2008h . 



Merloni et alJ ( 2004f ) expect a weak evolution of 
M BH /M sph cx (1+z) - 4 - - 6 , when fitting the to- 
tal stellar mass and star formation rate density 
as a function of redshift and comparing that to 
the hard X-ray selected quasar luminosity func- 
tion, assuming that BHs only grow through accre- 
tion. Such a s l ope i s in agreement with work by 
Hopkins et al. ( 2009h who combine prior observa- 



tional constraints in halo occupation models with 
libraries of high-resolution hydrodynamic simu- 
lations of galax y merg ers. Using semi-analytic 
models, ICrotonl ( 2006f ) predicts an evolution of 
M B n/M sph cx (l+z) ° 4 ~ 12 . A more r apid e volu- 
tion is predicted by IWvithe fc Loebl |2003) who 
assume a self-regulated BH growth model and find 
Mbh/ M sp h oc (1+z) 15 , similar to our observa- 
tional result. 

However, the great advantage of the study 
presented here are the high-quality images at 
hand, allowing for a detailed bulge-to-disk de- 
composition of the host galaxy of the low- and 
intermediate-z Seyfert-1 galaxies. Combining data 
from a large sample of active galaxies, covering 
a redshift range from the local Universe out to 
z=4.5, all treated in a consistent manner, re- 
sults in smaller error bars on the predicted evo- 
lution than previous studies. Moreover, it al- 
lows, for the first time, to correct evolutionary 
trends for selection effects. The evolution we find 
(M BH /M sph cx(l + z) 14±0 - 2 ) is indicative of BH 
growth preceding host-spheroid assembly. 

Still, we did not take into account that the evo- 
lution may depend on BH mass (see also Paper 
III). Indeed, there are theoretical predictions that 
objects with higher BH (or bulg e) masses evolve 



faster (e.g.. iHopkins et al.M2009l ). For example, 
iDi Matteo et all (|2008h find that when restricting 
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their fits to objects with M* > 5 X 10 10 M , the 
relation has a slope of ~ 1.9 at z=3-4 and ~ 1.5 
at z=2. Unfortunately, our sample is too small to 
allow us to address this possibility. There may in- 
deed be some evidence that the offset in BH mass 
is larger for objects with more massive BHs (Fig.[5j 
upper left panel and Fig. right panel). 

8.3. BH Mass - Galaxy Luminosity Rela- 
tion 

A different scenario seems to emerge when con- 
sidering the relation between Mbh and total host- 
galaxy luminosity (Fig. [5J lower left panel). This 
relation is almost non-evolvin g within the last 
six billion years Recently, Ijahnke et al. (2009) 
found qualitatively similar results for a small sam- 
ple of ten AGNs at redshifts between 1 < z < 2: 
They derive host-galaxy masses from colors based 
on ACS and NICMOS imaging, finding that they 
lie on the MTm-M+.h.,,t„p. re lation in the local Uni- 
verse (|Haring k RiJ2004 l. 

Such a non-evolving MsH-^host relation can be 
interpreted twofold. 

(a) The amount by which some of the more 
distant objects have to grow their spheroid is 
already contained within the galaxy itself, and 
the growth can be achieved by the redistribution 
of stars, i.e. transforming disk stars into bulge 
stars. Such a redistribution can be the result of 
mergers or secular evolution, e.g. bar instabili- 



ties (e. 



1998: A vila- Reese et al.l 



Combes k SandersHl98lHvan den Bosch 
20051 : iDebattista et al 



20061 ) an d torque-driven accretion (see e.g. re- 
view in iKormendv k KennicuttJ 2004 ) which 
may coincidentally be also the triggering mech- 
anism for the BH activit y we observe (e.g., 
Shlosman k Noeuchi 11993c lAthanassoula et al 



2003tlDumas et alJl2007tlHaan et alJl2009f ). While 
not every object in the intermediate-z sample will 
experience a major merger in the last 4-6 billion 
years, secular evolution is a promising alternative 
way to grow the spheroidal components in these 
objects. But even if they do experience a ma- 
jor merger (as indeed evidenced for at least some 
objects in our sample), the role of the merger 
depends on the merger type as discussed above 
(e.g. a merger between similar objects will simply 



move the system along the local relation). There 
may again be a dependency on BH mass: For 
the low-mass objects, the offset becomes almost 
negative, indicating that in the low-mass range, 
either the BH is, at the same time, still growing 
by a non-negligible amount (consistent with the 
higher Eddington ratio in the low-mass regimJ^l) 
or that not all of the stellar mass will end up in 
the spheroid component. Indeed, for local RM 
AGNs, at least 6/19 objects reside in late-type 
host galaxies (preferentially those with lower BH 
masses) . 

(b) The relation between BH mass and host- 
galaxy luminosity (or mass) may be the mor e fun- 
damen tal one. Indeed, this is predicted by Peng 
(|20071 ): In his thought experiment, he shows that a 
tight linear relation between Mbh and host-galaxy 
mass can evolve - if the galaxy mass function de- 
clines with increasing mass - due to "a central- 
limit-likc tendency for galaxy mergers, which is 
much stronger for major mergers than for minor 
mergers, and a convergence toward a linear rela- 
tion that is due mainly to minor mergers" . Also, it 
is possible that BHs in late-type galaxies or galax- 
ies without classical bulges, while not following 
the same Mbh scaling relations as spheroids (see 
discussion in § I8.2.1|) . they instead obey a more 
fundamental relation between BH mass and host- 
galaxy mass. 

However, the relation between host-galaxy lu- 
minosity and Mbh seems to exist only up to z<l: 
The offset for the high-z comparison s ample does 



not dec rease as the luminosity given bv lPeng et al 



(2006b) is already the total host-galaxy lumi- 
nosity (the same is true also for the r esults by 
Merloni et all 120091: iDecarli et all l2009t) . Along 



the line of argument of (a) above, the growth 
of the spheroid above a redshift of z > 1 can- 
not simply be achieved through secular evolution 
(with quasars being predominantly hosted by el- 
lipticals), but instead, major mergers are needed. 
A major merger is more likely to happen for the 
high-z sample given the longer time span. Or, fol- 
lowing (b), a relation between BH mass and host 
galaxy is already at place at z < 1 , but still evolv- 
ing at earlier times. However, we cannot exclude 
that part of the difference is due to the difference 



12 Note that there is insufficient information to constrain the 
intrinsic scatter. 



13 Of course, we may be biased against low-mass objects with 
low Eddington ratios. 
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iii BH mass between the samples, with the high-z 
objects generally having larger BH masses. 

In the end, the discussion boils down to the 
following question: What is the dominant mech- 
anism that grows spheroids, and does it depend 
on spheroid mass and/or redshift? This is de- 
bated controversially in the theoretical literature. 
For example , base d on their semi-analytic models, 
Parry et al.l ( 2009t ) find that the majority of ellip- 
ticals and spirals never experience a major merger 
but rather, that they acquire their spheroid stel- 
lar m ass through minor me rgers or disc instabil- 
ities. Hopkins et al.1 ( 2009c ). on the other hand, 
combine empirically constrained halo occupation 
distributions with high-resolution merger simula- 
tions, and find that major mergers dominate the 
formation of ^L* bulges and systems with higher 
B/T, but that lower- mass or lower B/T systems 
are preferentially formed by minor mergers. They 
predict that the major merger rate increases with 
redshift. Qualitatively, we can reconcile such a 
scenario with our results: Higher-mass objects and 
those at higher rcdshifts (i.e. the majority of the 
high-z sample) form their spheroids preferentially 
through major mergers and are thus still evolv- 
ing toward a MeH-ihost relation, while lower-mass 
and lower-z objects (i.e. our intermediate-z sam- 
ple) grow their spheroids through minor mergers 
or disk instabilities that redistribute the stars and 
thus, they fall on the AfeH-ihost relation. 

9. SUMMARY 

We study the evolution and intrinsic scatter of 
the BH mass - spheroid luminosity relation, taking 
into account selection effects, by combining three 
different samples of AGNs. Our intermediate- 
redshift sample comprises 40 Seyfert-1 galaxies at 
two different redshift bins (34 objects at z=0.36, 
and 6 objects at z=0.57; look-back time 4-6 Gyr) 
for which we measure the BH mass from single- 
epoch Keck spectra. The sample spans more 
than one order of magnitude in BH mass (log 
M BH /M =7.5-8.8). 2D surface-brightness pho- 
tometry using GALFIT is carried out on high- 
resolution HST images to decompose the image 
into AGN and host-galaxy components. The low- 
z comparison sample consists of 19 local AGNs 
(0.02 < z < 0.29; z avs = 0.08) with reverber- 
ation BH masses ( Bentz et al. 2009bl ). We re- 



analyzed the archival HST images in a way com- 
parable to our intermediate-z Seyfert galaxies to 
eliminate possible systematic offsets. Finally, we 
combine our results with high-z data (44 quasars 
from 0.66 < z < 4.5; z avo = 1.8) compiled from 
the literature, mainly consisting of gravitationally- 
lensed AGNs (jPeng et alJ2006bh that were treated 
in a self-consistent manner. For all objects, the 
spheroid luminosity is corrected for passive lumi- 
nosity evolution. Our main results can be summa- 
rized as follows. 

• We determine the evolution in Mbh with an 
unprecedented accuracy, taking into account 
selection effects. Our result, MBH/i sp hOc(l+ 
z) 1 indicates that BH growth precedes 
host-spheroid assembly. The intrinsic scat- 
ter, assumed to be non-evolving, is non- 
negligible (0.3±0.1 dex; <0.6 dex at 95% 
CL). It may reflect the different ways and 
timescales involved when growing spheroids 
or may partially be due to a high fraction of 
spirals and/or potential pseudobulges in our 
sample. 

• The local and intermediate-z sample fol- 
low an apparently non-evolving Mbh- host- 
galaxy luminosity relation. Either the 
spheroid grows by a redistribution of stars, 
or the relation between BH mass and host 
galaxy is more fundamental. Above z ~ 
1, the relation seems to be still forming, 
e.g. through major mergers. 

We are currently studying the evolution of the 
BH mass - spheroid velocity dispersion relation 
(Woo et al. 2009, in preparation), which should 
allow us to tighten the error bars on evolution 
given that velocity dispersion can be measured 
more precisely than host luminosity. Studying 
this independent relation will also enable us to 
distinguish between different evolutionary scenar- 
ios, probe the "fun damental plane" betw een Mbh, 
i sp h, and a (e.g., Hopkins et al. 20071 ). and per- 
form further tests for systematics. Due to the 
failure of NICMOS in Fall 2008, nine objects at 
z = 0.57 and three objects at z = 0.36 were 
not observed. Instead, we were recently allocated 
time with WFC3 to complete the full sample of 
Seyfert-1 galaxies. We will present results for this 
extended sample in another paper. At the same 
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time, increasing the local AGN comparison sample 
would be desirable (and indeed an HST proposal 
for the eight nearest RM AGNs by Bentz et al. is 
in the queue). Understanding slope and scatter of 
the local relations for active galaxies is crucial to 
study their evolution. 
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A. MONTE CARLO SIMULATIONS 



To probe the reliability of GALFIT to derive the AGN and host-galaxy properties accurately and to 
estimate the systematic uncertainties involved in the fitting, we ra n Monte C arlo simulations of a set of 
different galaxy models. A comparable procedure was carried out bv lKim et al. I (|2008l) . 

In particular, we used GALFIT to simulate galaxies, consisting of (a) PSF plus spheroid, and (b) PSF 
plus spheroid plus disk using a range of typical galaxy properties of our sample. In both cases, the total 
magnitude was set to either 18 or 19 mag. For case (a), we assumed the effective radius of the spheroid to 
be r e ff = 4 pix, 6 pix, 8 pix, and an AGN-to-total luminosity ratio of 0.05, 0.1, 0.2, 0.5, 0.8, 0.9, 0.95. In 
case (b) , the effective radius of the spheroid was set to r c g = 3 pix, 4 pix, 6 pix, 8 pix (to additionally probe 
the lower limit which can be a problem when fitting spheroid plus disk). For an AGN-to-total luminosity 
ratio of 0.8 (0.5), the spheroid-to-disk ratio was 0.5 (0.2, 0.5), and for an AGN-to-total luminosity ratio of 
0.2 and 0.1, a spheroid-to-disk ratio of 0.1, 0.2, 0.5 was used. 

Note that this is a conservative approach, focusing on the parameter space for which the detection of the 
spheroid component is most difficult, i.e. a small spheroid size, a small spheroid-to-disk luminosity ratio and 
a large AGN-total galaxy luminosity ratio. The galaxies were simulated with a given PSF and noise was 
added based on a typical observed signal-to-noise (SN) ratio in a Monte Carlo fashion creating 100 artificial 
images which were then fitted by GALFIT. In a first run, the simulated galaxies were fitted with the same 
PSF that was used to create the artificial images, in subsequent runs with different PSFs from our PSF 
library to simulate PSF mismatch. 

For case (a), GALFIT can easily recover the sizes and magnitudes, even when the spheroid reaches sizes 
close to the minimum size that can be resolved given the PSF (here assumed to be 2.5 pixels). In Fig. [8j 
we show the resulting offsets for the smallest spheroid (r e g = 4 pix). However, more caution needs to be 
exercised for a three-component fit (PSF, spheroid, and disk), probed in case (b). In Fig. [9) we show the 
resulting offsets for the smallest spheroid (r c g — 3 pix), i.e. the most difficult scenario for retrieving the 
spheroid parameters accurately. The derived spheroid magnitude can differ up to 0.5 mag in the worst case, 
while the difference in PSF magnitude is less than 0.2 mag. We adopt these values as conservative measures 
of our errors, i.e. 0.2 mag for AGN luminosity and 0.5 mag for spheroid luminosity. 

As the estimation of errors is the main purpose of thi s analysis, we do not further discuss the results 



of these simulations. The overall trend is the same as in iKim et alj (|2008l ). i.e. the scatter in all derived 



parameters is largest when the AGN is dominant, and when r e g is small and difficult to distinguish from 
the nucleus or large with low surface brightness. Spheroid-to-disk-to-AGN decompositions are much more 
difficult than spheroid -to-AGN as they involve 6 additional free parameters (if the spheroid is fitted by a 



dimcuit tnan spncroid-to-AUiN as tney involve o additional tree paran 
de Vaucouleur s (1948) profile) and can only be done if the S/N is high. 



B. CHOICE OF SERSIC INDEX 

A general profile to fit galaxies is the so-called [Sersic (1968) power law, which is defined as 



E(r) = S off exp 




where £ e ff is the pixel surface brightness at the effective radius r e g 



(Bl) 

and n is the Sersic index. In this 

generalized form, an exponential disk profile has n = 1, a de VaucouleursI ( 19481 ) p rofile has n = 4 , and a 
Gaussian has n = 0.5 (which was used in Paper II to fit a bar component). While de VaucouleursI ( 19481 ) 
profiles are traditionally and widely used to fit spheroidal components, recent studies show that spheroids can 
have Sersic indices ranging between 0.5 and 6. Disk galaxies typica lly have a bulge comp onent with n < 4, 
with classical bulges having n > 2 and pseudobulges having n < 2 ( Fisher fc Drorv 2008). Moreover, there 
seems to be a relation betw e en the Sersic index and the spheroid luminosity or host-galaxy luminosity (e.g., 
Kormendv fc Bruzuall fl978l : IShaw fc Girmorelll989t lAndredakis fc Sandersj[l994t lGrahamll200ll iPeng et al 
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20021: iNipoti et al.|[20o5 iMacArthur et al.|[2008h . 

The objects we are fitting are complex in nature, in particular due to the presence of the AGN, a very 
luminous point source in the center for which a perfectly matching PSF fit cannot always be achieved. Thus, 
we cannot use a fit with the Sersic index n as a free parameter, as it would add yet another free component to 
an already difficult fit and increase degeneracies between PSF, bulge, and disk. Such an approach could easily 
lead to an unphysical fit, if GALFIT is trying to fit any remaining PSF mismatch with such a component. 
In such a situation, an alternative approach to estimate the best fitting Sersic exponent is to use a range of 
Sersic indices, keep them fixed at each step and then obtain the best n from the resulting the \ 2 statistics. 
This approach is generally recommen ded when attemp t ing galaxy deco mpositions of faint or difficult to 
model galaxies like AGN host galaxies ( Peng et al. 20021 : Kim et al. 2008 ). 

To test the systematic uncertainties in derived spheroid and PSF magnitude depending on the adopted 
Sersic index, we re-ran our models using n — 0.5, 1,2,3,4,5 to fit the spheroid component. At each step, 
we kept n fixed to the chosen value but allowed all other parameters to vary, including the disk component. 
For those objects that were initially (i.e. when using n=A) fit by a spheroid component only, we carefully 
checked the residuals of the resulting best n fit for any evidence of an additional disk component. For only 
one source (1501) is the quality of the fit increased significantly by the addition of a disk component 1 14 l 
Note that using the best-fitting n instead of n=4 docs not in general solve the problem of a a vanishingly 
small bulge component for some objects. For 11 of the 16 objects discussed in § I4.3[ nothing changes. For 
five objects, the effective radius of the bulge component is no longer smaller than the FWHM of the PSF; 
however, for two different objects, the bulge component then becomes vanishingly small. 

The results are shown in Fig. [Hjl separating objects for which the host galaxy was fitted by a spheroid 
component only (left panel) and those for which the host was fitted by a spheroid plus disk component (right 
panel). While the overall trend is the same, fitting the host galaxy by two components results in a larger 
scatter because the disk magnitude can also vary. The results can be summarized as follows: Decreasing n 
from 4 to 3 (2, 1, 0.5) decreases the spheroid luminosity - on average by 0.08 (0.23, 0.41, 0.54) mag - and 
increases the nuclear luminosity - on average by -0.07 (-0.24, -0.37, -0.4) mag. Increasing n from 4 to 5 (6) on 
the other hand increases the spheroid luminosity on average by -0.07 (-0.13) mag and decreases the nuclear 
luminosity on average by 0.1 (0.22) mag. For all but the most extreme indices, potential systematics related 
to the choice of Sersic index are small compared to the adopted uncertainty on the spheroid luminosity (0.5 
mag) and on Mbh (0.4 dex). 

Another approach is to calcul ate the best Sersic i ndex we would expect based on the measured host-galaxy 
magnitude using the relation in Nipoti et alJ ( 2008 ), de r ived from surface -photometry of ACS images of the 



well-defined Virgo cluster sample ( Ferrarese et al. 2006 : Gallo et al. 20081 ) 



logn = (0.27 ±0.02) log 



L 



sph,B 



9.27 +0.4 ±0.02 



(B2) 



For the 23 objects studied here, we estimate a Sersic index ranging from ~4.0 to 5.9, on average 4.9 ± 0.5. 
For the sample studied in Paper II, the Sersic index ranges from ^4 to 6.8, with an average of 5.5 ± 0.6. 
Within the errors, these values are in agreement with estimates using the relation between Sersic index and 
bulge B-band magnitude for a local sample from iGrahaml (|200lh (their Figure 14, middle panel). Note that 
in both cases, the estimated Se rsic index remains the same within the errors when using t he host-galax y 
luminosity (for n-L relation from lNipoti et al1l2008h or bulge luminosity (for n-L relation from lGrahaml l200£ 
as derived from the best fit with a free Sersic index instead of the one derived from n=4. As the relation 
between Sersic index and host-galaxy magnitude has its own uncertainties and scatter, and as our average 
value is close to 4, we adopted the simpler solution of fixing n to 4 for all objects as our default choice. This 
also allows a better comparison with other AGN host-galaxy studies. 

To ultimately probe the potential systematics related to the choice of Sersic index, we performed the 
same analysis as for n=4, but this time using the best n derived from the procedure described above (i.e. as 



*Note, however, that fitting this source is particularly complicated as it is in the process of merging with a neighboring galaxy. 
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chosen based on the \ 2 statistics when performing a variety of fits with n fixed to 0.5, 1, 2, 3, 4, 5, 6) both 
for the intermediate-z and the local sample. The same approach was followed for the local sample of RM 
AGNs. None of the results stated in the paper change: The resulting fits, offsets, and predicted evolutionary 
trends remain the same within the errors. More precisely, for the evolution in AIbh(Mbh/ L sp hOc(l + z) 13 ), 
including selection effects, we obtain = 1.3±0.2 (instead of = 1.4±0.2 for n — 4) for the full sample with 
an intrinsic scatter <0.7 dex at 95% CL (0.4±0.1 dex for a uniform prior on (Tint) and = 3. 4 ±1.2 (instead 
of = 2.8 ± 1.2 for n = 4) for the intermediate sample alone, adopting again the prior by IGiiltekin et all 
(|2009h . 



SURFACE PHOTOMETRY OF RM AGNs 



For an homogen eous treatment of al l data, we performed an independent analysis of the HST archival 
images presented in lBentz et al. I (l2009bh . u sing the sam e appro ach as for our distant Seyfert galaxies (§EJ). 
Details of the observations can be found in Bentz et alj ( 2009bl ). 

We disregarded the five objects observed with WFPC2/PC due to the low quality of the data, the PSF 
mismatch when using a synthetic PSF created by TinyTim and the lack of stellar PSFs on the images. These 
problems made it difficult to achieve satisfactory fits. The spheroid radius found with GALFIT was either 
in the lower limit of 2.5 pixels (=FWHM of PSF) - probably because it was fitting a PSF mismatch - or 
was unphysically large. Thus, we here focus on the ACS/HRC data alone. 

From the remaining 30 objects imaged with ACS/HRC, we first excluded all NGC objects (8/30) which 
are nearby and extended and for which the field-of-view is too small to measure the sky background. Also, 
they are often affected by dust lanes. The latter is also the case for IC4239A (plu s an unreliab l e BH m ass). 
For the same reasons, these objects were also excluded in the further analysis bv lBentz et al. ( 200 9a). We 
decided to additionally exclude Fairall 9 due to a dust lane crossing the spheroid and PG0953+414 for which 
no reasonable fit could be achieved. Thus, our final robust sample consists of 19 objects. 

We used the pipeline-processed data and combined them using multidrizzle, to remove cosmic rays and 
defects and correct for distortion. (Note that multidrizzle takes into account the saturated pixels of the 
longer exposures and combines the images accordingly.) As the data were not dithered, no improvement 
of sampling was achieved and the final scale is 0.025 arcsec/pixel (pixfrac=0.9). For these ACS/HRC data 
imaged in the F550M filter, the PSF created by TinyTim is not as good a match as it is for the NICMOS 
images. We therefore additionally created a PSF from a star observed in one of the images (Mrk 110) and 
performed extensive tests to compare their quality. As the TinyTim PSF typically gave a bad fit in the core, 
but the stellar PSF had too low S/N in the wings, we decided to combine both PSFs (the synthetic PSF 
for the wings, the stellar PSF in center out to r — 2 x FWHM), which significantly improved the quality 
of the fits. This PSF enabled us to fit the AGN with on ly one PSF without the need of corrections of PSF 



mismatch (e.g. by using an additional PSF as done by IBentz et al. lj2009tjn . We used the same criteria 
as for our distant Seyfert sample to d ecide whether we need to fit an additional disk component (see § |4j . 



For four objects, IBentz et al. ( 2009b ) fitted both a spheroidal and disk component, while we decided that 
fitting a spheroidal component alone is sufficient. One object has a saturated PSF (PG1226+023) and we 
masked out the saturated center to fit the PSF to the wings only. For three objects (Ark 120, Mrk 279, 
and PG 1211+143), we out-masked the very center of the PSF and fitted the PSF to the wings only due to 
remaining PSF mismatch. 

We compare the results in Fig. [TT1 For this compa rison, we add the different PSF components and the 



different spheroid components of IBentz et al 



( 200911 to a "total" PSF magnitude and "total" spheroid 
magnitude, respectively. While the PSF and tota l magnit udes generally agree well, the spheroid magnitudes 
we derive are often brighter than in IBentz et al Ij20p9b), especia lly in those cases where we fit a spheroid 
component only and not spheroid+disk as in IBentz et al.l ( |2009bl ) (4 objects). 

As for our intermediate redshift sample (see Appendix [BJ, we also calculated the best Ser s ic ind ex we 
would expect based on the measured host-galaxy magnitude using the relation in lNipoti et al. I (|2008i ). The 
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average value of n=5.2 ± 1.6 is in agreement within the errors wit h the ave r age n derived when using 
the relation between Sersic index and bulge B-band magnitude from iGrahaml ( 2001 ). It also agrees well 
with the average n estimated for the intermediate-z sample, with a larger scatter due to the larger spread 
in luminosities. We carefully checked whether when using the best-fitting Scrsic index, there is the need of 
adding a disk component for those objects for which the host galaxy was originally fitted by a n—4 component 
only; we do not find such evidence in any of the objects. 
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Fig. 1. — Postage stamp NICMOS images of the 23 ^cyfert-1 galaxies in the sample. The first 17 objects 
are at z = 0.36, the last six objects are at z = 0.57. A 4 arcsecond scalebar is shown in the upper left image, 
corresponding to <~ 20 kpc at z=0.36 and 26 kpc at z=0.57, respectively. The label M or M/l marks objects 
that are apparently merging or interacting. 




semi-major axis (arcsec) 



Fig. 2. — Surface-brightness profiles for all objects with a resolved spheroid component, measured from the 
data as well as from each component that was fitted. Note that the fits were performed in two dimensions 
using GALFIT, so this figure is for illustration purposes only, showing the relative contribution of each 
component as a function of radius. Some profiles show an early truncation which is an artifact of the 
elliptical isophote routine used to make the plots due to nearby objects. (For the measurements, these 
objects were fitted simultaneously using GALFIT.) [See the electronic edition of the Journal for a color 
version of this figure.] 
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Fig. 3. — The same as in Figure [2] for objects for which the bulge models correspond to the minimum size 
allowed by HST resolution. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 4. — Histogram of the bulge-to-total flux ratios (B/T) from iBenson et al.l ([20071 ) for galaxies within 
0.5 mag of the Seyfert total host-galaxy magnitude (black line). The red dashed line shows the upper limit 
on B/T we derived for the Seyfert galaxies from GALFIT. We use this upper limit to cut the distribution 
and to calculate a mean (blue line) and sigma of the remaining B/T values (red shaded area). The first 11 
objects were imaged with NICMOS (ten at z—0.36, one at z=0.57) and the last 5 objects were studied in 
Paper II, but we include them here to estimate spheroid luminosities from upper limits. [See the electronic 
edition of the Journal for a color version of this figure.] 
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Fig. 5. — Upper left panel: Black hole mass-spheroid V-band luminosity relation. Colored circles represent 
measurements for the intermediate-redshift Seyfert galaxies (red: z = 0.57, green: z = 0.36, blue: z = 0.36 
taken from Paper II; squares indicate objects for which the fitting procedure ran into the lower limit of 
the spheroid effective radius and we used priors to obtain a meas ure of the spheroid luminosity). Black 
circles correspond to the local RM sample (z avc ~ 0.08) studied bv lBentz et al. ( 2009a. b]) and re-analyzed 
here, including the best fit (black solid line; see text and Table Q] for details). For all objects, the spheroid 
luminosity is evolved to z = assuming pure luminosity evolution (see text for details) . Note that no selection 
effects are included here. Intermediate-z objects with signatures of interaction or mergers (see Fig. [T] and 
Paper II) are in dicated by a large open black circle. The dashed line shows the fiducial local relation for 
inactive galaxies ( Marconi fc Hunt! 20031 ). transformed to V-band (group 1 only; see text for details). Upper 
middle panel: The same as in the left panel, th is time all z = 0.36 objects in blue. Green circles are the 
high-z AGN sample (average z ~ 1.8) taken from lPeng et al. ( 2006bh and treated in a comparable manner. 
We assume 0.4 dex as error on -Mbh, and 0.12 dex as error on luminosity (based on the error quoted by 
Peng et all (|2006bh of 0.3 mag). We mark those high-z objects for which the BH mass is based on the C IV 
line as green squares. Upper right panel: Distribution of residuals in log Mbh with respect to the fiducial 
local relation of RM AGNs. Top panel: distribution of residua ls for intermediate- redshift Seyfert galaxies 
(blue: z=0.36; red: z=0.57) and for the high-z AGN sample from Peng et al. ( 2006bf ) (green). Bottom panel: 
local sample. Lower panels: The same as in the upper panels, for the total host-galaxy luminosity. [See 
the electronic edition of the Journal for a color version of this figure.] 
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Fig. 6. — Left panel: Offset in log Mbh as a function of log (1 + z) with respect to the fiducial local relation 
of RM AGNs (Fig. El upper middle panel). The best fit to all data points (solid black line) of the form 
A log Mbh = jlog(l + z) including intrinsic scatter in log Mbh as a free parameter but ignoring selection 
effects is 7 = 1.2 ± 0.2. (Note that the average data points for each sample are plotted only to guide the 
eye.) For comparison, we also overplot the selection-bias corrected evolution (Mbh/AspIi oc (1 + z) 1 °- 2 ; 
dotted line) with the la range as dashed lines. As in Fig. [5l squares indicate objects for which the fitting 
procedure ran into the lower limit of the spheroid effective radius and we used priors to obtain a measure of 
the spheroid luminosity. Right panel: The same as in the left panel as a function of look-back time. Here, 
the symbol size corresponds to BH mass. 
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Fig. 7. — Results of Monte Carlo simulations probing the effect of selection effects on the slope (3 of the 
relation A log Mbh= /31og(l 4- z) at fixed zero redshift spheroid luminosity corrected for evolution, and 
intrinsic scatter a- m t of the MeH-isph relation which is assumed to be non-evolving. Plotted are the 68% and 
95% joint confidence contours. Left panel: Including both intermediate-z and high-z sample, without an 
assumed prior on <7i nt . Both and a- m t are well const rained (0 = 1.4 ± 0.2; gj nt = 0.3 ± 0.1). Middle panel: 
The same as in the left panel, including the prior bv lGultekin et al~ ( 20091 ) (i.e. <7i n t = 0.38±0.09), resulting 
in the same within the errors. Right panel: The same as in the middle panel, but for intermediate-z 
sample only. Whi le our sample a l one d oes not cover a large enough range in redshift, we find /3=2.8±1.2 
using the prior bv lGultekin et al. ( 20091 ) on a mt . 
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Fig. 8. — Results of GALFIT fits to simulated galaxies, consisting of PSF plus spheroid, with noise added in a 
Monte Carlo fashion, realizing 100 artificial images for each parameter combination. The difference between 
input AGN magnitude and derived AGN magnitude is shown (upper panels), the difference between input 
spheroid magnitude and derived spheroid magnitude (middle panel), and the difference between input and 
derived effective radius of the spheroid (lower panel). Each data point represents the average plus error of 
GALFIT fits to the 100 artificial images. Black data points correspond to fits where the PSF used to create 
the artificial image is identical with the one used for fitting. The other three colors correspond to a different 
PSF used for fitting which was taken from our PSF library to simulate PSF mismatch. The left figure shows 
results for a total host-galaxy magnitude of 18 mag, the right figure of 19 mag. The left panels within each 
figure correspond to an axis ratio of b/a = 0.5, the right panels to b/a = 0.1. In both figures, the input 
effective radius of the spheroid component is set to 4 pixels. 
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Fig. 9. — The same as in Fig. [8] for artificial images consisting of PSF plus spheroid plus exponential 
component. The difference between input disk magnitude and derived disk magnitude and the difference 
between input disk radius and derived disk radius is shown additionally in the two lower panels. In the left 
figure, the spheroid-to-disk ratio is 0.2 (and thus, the AGN-to-total luminosity /agn plotted on the x-axis 
only assumes values of 0.1, 0.2, 0.5), in the right figure, the spheroid-to-disk ratio is 0.5 (with /agn = 0.1, 
0.2, 0.5, 0.8; see text for details). In both cases, the total host-galaxy magnitude is 18 mag and the effective 
radius of the spheroid is set to 3 pixels. While the PSF magnitude can be retrieved easily to within 0.2 mag, 
the difference in spheroid magnitude can be up to 0.5 mag in the worst case. 



34 




123456 123456 
sersic index sersic index 



Fig. 10. — Systematic effects in derived magnitudes of spheroidal and AGN component due to the adopted 
spheroid profile (Sersic index 0.5, 1, 2, 3, 4, 5, or 6). Left panel: All objects for which the host galaxy was 
fitted by a spheroid component only, with each object corresponding to a given color. Right panel: Same 
as the in the left, for objects for which the host galaxy was fitted with a spheroid plus disk component. See 
text for details. [See the electronic edition of the Journal for a color version of this figure.} 
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Fig. 11. — Difference between the results of the surface-brightness fitting of iBentz et al. (2009b) and our 
work here, for PSF magnitude (left panel), spheroid magnitude (middle panel), and total magnitude 
(right panel). 
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Table 1 
Sample Properties 



Name 




z 


Dl 


RA (J2000) 


DEC (J2000) 


i 










Mpc 










mag 


(i) 




(2) 


(3) 




(4) 




(5) 


(6) 


0107 


(Sll) 





3558 


1892.9 


01 


07 15 


97 


-08 34 29.4 


18.47 


0804 


(SSI) 





3566 


1897.9 


08 


04 27 


99 


+52 23 06.2 


18.55 


0934 


(SS2) 





3672 


1964.1 


09 


34 55 


60 


+05 14 09.1 


18.82 


1007 


(SS5) 





3733 


2002.5 


10 


07 06 


26 


+08 42 28.4 


18.69 


1015 


(S31) 





3505 


1860.0 


10 


15 27 


26 


+62 59 11.5 


18.15 


1021 


(SS6) 





3584 


1909.1 


10 


21 03 


58 


+30 47 55.9 


18.92 


1043 


(SS7) 





3618 


1930.3 


10 


43 31 


50 


-01 07 32.8 


18.82 


1046 


(SS8) 





3656 


1954.1 


10 


46 10 


60 


+03 50 31.2 


18.45 


1258 


(SS9) 





3701 


1982.3 


12 


58 38 


71 


+45 55 15.5 


18.56 


1334 


(SS10) 





3658 


1955.4 


13 


34 14 


84 


+ 11 42 21.5 


17.83 


1352 


(SS11) 





3732 


2001.8 


13 


52 26 


90 


+39 24 26.8 


18.39 


1501 


(SS12) 





3625 


1934.7 


15 


01 16 


83 


+53 31 02.4 


17.80 


1505 


(SS13) 





3745 


2010.0 


15 


05 41 


79 


+49 35 20.0 


18.73 


1611 


(S28) 





3679 


1968.5 


16 


11 56 


30 


+45 16 11.0 


18.63 


2115 


(SS14) 





3706 


1985.5 


21 


15 31 


68 


-07 26 27.5 


19.24 


2158 


(S29) 





3575 


1903.5 


21 


58 41 


93 


-01 15 00.3 


18.95 


2340 


(SS18) 





3582 


1907.9 


23 


40 50 


52 


+01 06 35.5 


18.50 


0155 


(Wll) 


0.5634 


3270.9 


01 55 16.18 


-09 45 56.0 


20.09 


0342 


(W22) 


0.5648 


3280.8 


03 42 29.70 


-05 23 19.5 


18.70 


1439 


(W12) 


0.5623 


3263.2 


14 39 55.11 


+35 53 05.4 


19.02 


1500 


(W20) 


0.5753 


3354.7 


15 00 14.81 


+32 29 40.4 


19.60 


1526 


(W16) 


0.5782 


3375.2 


15 26 54.93 


-00 32 43.3 


19.99 


1632 


(W8) 


0.5703 


3319.4 


16 32 52.42 


+26 37 49.1 


18.70 



Note. — Col. (1): Target ID (RA: hhmm). In brackets, the name used in 
other publications. Col. (2): Rcdshift from SDSS-DR7. Col. (3): Luminos- 
ity distance in Mpc, based on rcdshift and the adapted cosmology. Col. (4): 
Right Ascension. Col. (5): Declination. Col. (6): Extinction-corrected i' 
AB magnitude from SDSS-DR7 photometry ( "modclMag_i" ) . 
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Table 2 

Results from Imaging of Distant Seyfert Sample 



-a 



Name 


Total 


Host 


Spheroid 


lug ^host,V/ ^0 


lr>0" T, i t t I J / r-. 

lu & ^sph, V / -^O 
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■^-^5100 


J nuc 
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mag 
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(A\ 
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\ 1 ) 


[Pi 


(■') 


(HI) 


( LL > 


0059+1538 (S09; M) 


18.24 


18.51 


19.08+0.50 


10.91 


10.68 


3.24 


0.71 


0.22 


8.13 


3 


0101-0945 (S10) 


18.03 


18.37 


19.32+0.58 


10.95 


10.28(10.57) 


0.49 


1.03 


0.27 


8.25 


3 


0213+1347 (S12) 


18.20 


18.56 


21.23+0.60 


10.90 


9.56(9.83) 


0.54 


0.97 


0.28 


8.67 


3 


1105+0312 (S21; M) 


17.49 


17.94 


18.99+0.58 


11.13 


10.39(10.71) 


0.51 


2.15 


0.34 


8.79 


3 


1119+0056 (S16) 


19.16 


19.87 


22.28+0.50 


10.41 


9.45 


0.76 


0.73 


0.48 


8.27 


3 


1400-0108 (S23) 


18.02 


18.39 


20.88+0.55 


10.94 


9.72(9.95) 


0.57 


1.11 


0.29 


8.70 


4 


1400+0047 (S24) 


18.09 


18.22 


18.61+0.50 


11.05 


10.89 


12.65 


0.44 


0.11 


8.33 


3 


1529+5928 (S26) 


18.88 


19.22 


20.07+0.50 


10.67 


10.33 


0.75 


0.52 


0.27 


8.02 


3 


1536+5414 (S27; M/I) 


18.53 


19.00 


19.48+0.50 


10.75 


10.56 


4.78 


0.95 


0.36 


8.10 


3 


1539+0323 (SOI; M) 


18.54 


18.91 


19.97+0.50 


10.77 


10.34 


5.30 


0.72 


0.29 


8.20 


4 


1611+5131 (S02) 


19.04 


19.32 


19.87+0.50 


10.58 


10.36 


2.63 


0.34 


0.22 


7.98 


3 


1732+6117 (S03; M) 


17.97 


18.50 


20.25+0.53 


10.92 


9.97(10.22) 


0.50 


1.64 


0.39 


8.28 


4 


2102-0646 (S04) 


18.12 


18.61 


20.18+0.50 


10.88 


10.25 


0.96 


1.33 


0.36 


8.44 


4 


2104-0712 (S05) 


18.00 


18.68 


20.51+0.50 


10.83 


10.10 


1.03 


1.85 


0.47 


8.74 


4 


2120-0641 (S06) 


18.48 


18.70 


20.62+0.50 


10.88 


10.11 


1.01 


0.51 


0.18 


8.16 


4 


2309+0000 (S07; M/I) 


17.82 


18.48 


20.39+0.50 


10.91 


10.15 


1.01 


2.10 


0.45 


8.53 


3 


2359-0936 (S08) 


18.33 


18.89 


21.77+0.50 


10.77 


9.62 


1.23 


1.22 


0.40 


8.10 


4 


0107-0834 (Sll) 


17.85 


18.01 


18.84+0.50 


10.84 


10.51 


0.59 


0.52 


0.14 


8.00 


3 


0804+5223 (SSI) 


17.89 


18.01 


19.34+0.58 


10.84 


10.04(10.31) 


0.47 


0.39 


0.11 


7.75 


3 


0934+0514 (SS2; M/I) 


18.38 


18.53 


18.53+0.50 


10.67 


10.67 


2.55 


0.33 


0.13 


7.72 


2 


1007+0842 (SS5) 


18.34 


18.80 


19.69+0.65 


10.58 


9.89(10.22) 


0.49 


0.93 


0.34 


7.66 


3 


1015+6259 (S31) 


17.83 


17.91 


18.67+0.50 


10.86 


10.56 


1.07 


0.29 


0.08 


7.94 


3 


1021+3047 (SS6; M?) 


18.85 


19.14 


20.29+0.70 


10.40 


9.60(9.94) 


0.48 


0.37 


0.24 


7.47 


3 


1043-0107 (SS7) 


18.31 


18.45 


19.31+0.63 


10.68 


10.02(10.34) 


0.51 


0.31 


0.12 


7.53 


3 


1046+0350 (SS8) 


17.89 


18.04 


19.67+0.55 


10.86 


9.95(10.20) 


0.48 


0.51 


0.13 


7.89 


3 


1258+4555 (SS9) 


18.04 


18.37 


18.37+0.50 


10.74 


10.74 


1.62 


0.93 


0.26 


8.05 


2 


1334+1142 (SS10) 


17.58 


18.19 


18.68+0.65 


10.80 


10.25(10.60) 


0.48 


2.26 


0.43 


7.94 


3 


1352+3924 (SS11) 


18.13 


18.31 


19.44+0.60 


10.77 


10.02(10.32) 


0.49 


0.51 


0.15 


8.11 


3 


1501+5331 (SS12; M) 


17.38 


18.19 


18.19+0.58 


10.79 


10.69(10.79) 


0.48 


3.24 


0.52 


8.15 


2 


1505+4935 (SS13) 


18.40 


18.92 


18.92+0.50 


10.53 


10.53 


1.09 


0.98 


0.38 


7.63 


2 


1611+4516 (S28) 


18.08 


18.11 


18.86+0.50 


10.84 


10.54 


0.94 


0.11 


0.03 


7.90 


3 


2115-0726 (SS14) 


18.97 


19.20 


19.20+0.50 


10.41 


10.41 


1.59 


0.29 


0.19 


7.64 


2 


2158-0115 (S29; M?) 


18.36 


18.48 


19.18+0.65 


10.66 


10.05(10.38) 


0.48 


0.25 


0.10 


7.94 


3 


2340+0105 (SS18) 


18.41 


18.79 


20.20+0.65 


10.53 


9.66(9.97) 


0.48 


0.70 


0.30 


7.51 


3 


0155-0945 (Wll; M/I?) 


19.64 


19.82 


19.82+0.50 


10.64 


10.64 


2.17 


0.31 


0.15 


8.00 


2 


0342-0523 (W22; M/I?) 


18.05 


18.53 


18.53+0.50 


11.16 


11.16 


7.34 


3.17 


0.36 


8.65 


2 


1439+3553 (W12) 


18.53 


18.96 


19.21+0.65 


10.98 


10.54(10.88) 


0.62 


1.87 


0.33 


8.72 


3 


1500+3229 (W20) 


19.00 


19.16 


19.16+0.50 


10.93 


10.93 


3.08 


0.51 


0.13 


8.52 


2 


1526-0032 (W16; M) 


19.33 


19.58 


19.58+0.50 


10.76 


10.76 


1.33 


0.60 


0.21 


7.59 


2 


1632+2637 (W8) 


18.48 


19.08 


19.08+0.50 


10.95 


10.95 


1.52 


2.59 


0.42 


8.73 


2 



Note. — Col. (1): Target ID (RA: hhmm). In brackets, the name used in other publications. Additionally, M/I marks objects that arc merging/interacting 
(see Fig. [T]and Paper II). All S* and SS* ob jects are at z ~ 0.36, all W* objects are at 2 - 0.57. The first 17 objects were observed with ACS/F775W 
and are already included in Trcu ct al. (2007), but arc listed here again due to a small error in extinction correction (<0.15 mag); also the luminosity in V 



was not included in Trcu ct al. (2007). For those objects with upper limits in Paper II, we here estimate the spheroid luminosity using priors ( i|4.3t . Col. 
(2): Total extinction-corrected F110W AB magnitude (SS* and W* objects) or F775W AB magnitude (S* objects); including PSF (with an uncertainty 
of 0.2 mag). Col. (3): Total host-galaxy extinction-corrected F110W (F775W) AB magnitude (with an uncertainty of 0.1 mag). Col. (4): Spheroid 
extinction-corrected F110W (F775W) AB magnitude (with an uncertainty of 0.5 mag). Col. (5): Logarithm of total host-galaxy luminosity in rest-frame 
V (solar units), not corrected for evolution. Col. (6): Logarithm of spheroid luminosity in rest-frame V (solar units), not corrected for evolution. For those 
objects, for which the fitting procedure ran into the lower limit of the spheroid effective radius, we give the corresponding value as upper limit in brackets. 
Col. (7): Spheroid effective radius (in kpc; semi-major axis). Col. (8): Nuclear rest-frame luminosity at 5100A (in 10 44 ergs -1 ) (uncertainty of 20%). Col. 
(9): Nuclear light fraction in F110W (F775W) (uncertainty of 20%). Col. (10): Logarithm of BH mass (solar units) (uncertainty of 0.4 dex). For those 
objects, for which the fitting procedure ran into the lower limit of the spheroid effective radius, we give the corresponding value as upper limit in brackets. 
Col. (11): Number of components fitted (2=PSF+spheroid; 3=PSF+spheroid+disk; 4=PSF+spheroid+disk+bar) . 



Table 3 

Results from Imaging of Local Comparison Sample 



Name 




z 
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log Mbh/A/q 
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Mpc 
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(1) 
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(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


3C120 





03301 


144.9 


13.49 


14.62 


14.62 


10.41 


10.41 


3.26 


7.74±0.21 


2 




Q 


05610 


250 5 


14 90 


15 76 
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10 45 


10 45 


2 48 


o Aa-i-n i n 

O.^UXU. J-U 


2 


Arkl20 





03271 


143.6 


13.29 


14.13 


15.29 


10.59 


10.13 


0.09 


8.18±0.06 


3 


Mrk79 





02219 


96.7 


14.46 


14.96 


15.94 


9.91 


9.52 


0.83 


7.72±0.12 


3 


MrkllO 





03529 


155.2 


15.42 


16.22 


17.50 


9.83 


9.32 


0.37 


7.40±0.11 


3 


Mrk279 


0. 


03045 


133.5 


14.00 


14.90 


16.04 


10.22 


9.77 


0.56 


7.54±0.11 


3 


Mrk335 





02579 


112.6 


14.18 


15.32 


16.29 


9.90 


9.51 


0.45 


7.15±0.11 


3 


Mrk590 





02639 


115.3 


14.21 


14.24 


15.42 


10.35 


9.88 


1.10 


7.68±0.07 


3 


Mrk817 





03146 


138.0 


14.33 


14.97 


17.42 


10.22 


9.24 


0.08 


7.69±0.07 


3 


PG0052+251 





15500 


739.2 


15.02 


16.25 


16.25 


11.24 


11.24 


16.76 


8.57±0.09 


2 


PG0804+761 





10000 


460.3 


13.89 


16.50 


16.50 


10.71 


10.71 


3.73 


8.84±0.05 


2 


PG0844+349 





06400 


287.4 


14.28 


16.10 


16.10 


10.44 


10.44 


3.87 


7.97±0.18 


2 


PG1211 + 143 





08090 


367.6 


14.43 


16.93 


16.93 


10.33 


10.33 


3.06 


8.16±0.13 


2 


PG1226+023 





15834 


757.7 


12.86 


15.55 


15.55 


11.54 


11.54 


4.42 


8.95±0.09 


2 


PG1229+204 





06301 


282.7 


15.38 


15.66 


16.65 


10.60 


10.20 


1.24 


7.86±0.21 


3 


PG1411+442 





08960 


409.5 


14.58 


16.80 


16.80 


10.48 


10.48 


9.52 


8.65±0.14 


2 


PG1613+658 





12900 


605.2 


14.48 


15.48 


15.48 


11.37 


11.37 


19.54 


8.45±0.20 


2 


PG1700+518 





29200 


1505.1 


14.87 


17.84 


17.84 


11.41 


11.41 


15.79 


8.89±0.10 


2 


PG2130+099 





06298 


282.6 


14.64 


16.37 


17.87 


10.32 


9.72 


4.15 


7.58±0.17 


3 



Note. — Results from imaging of local comparison RM AGN sample. Details of observations arc given in Bcntz ct al. (2009b). Briefly, all 
objects considered here were imaged with HST/ACS, in the F550M filter using the HRC chip. 

Col. (1): Target ID. Col. (2): Rcdshift. Col. (3) Luminosity distance in Mpc, based on rcdshift and the adapted cosmology. Col. (4): 
Total extinction-corrected F550M AB magnitude, including PSF (uncertainty of 0.2 mag). Col. (5): Total host-galaxy extinction-corrected 
F550M AB magnitude (uncertainty of 0.1 mag). Col. (6): Spheroid extinction-corrected F550M AB magnitude (uncertainty of 0.5 mag). Col. 
(7): Logarithm of total host-galaxy luminosity in rest-frame V (solar units), not corrected for evolution. Col. (8): Logarithm of spheroid 
luminosity in rest-frame V (solar units), not corrected for evolution. Col. (9): Spheroid effective radius (in kpc; semi-major axis). Col. (10): 
Logarithm of BH mass (solar units) with error, taken from lBentz et al.l l !2009bT) . Col. (11): Number of components fitted (2=PSF+sphcroid; 
3=PSF+spheroid+disk) . 



Table 4 

Fits to the local RM AGN log M B h - log £ sp h,v relation 



Method 


Sample 


K 


a 


Scatter 


(1) 


(2) 


(3) 


(4) 


(5) 


linear fit 


this work (with evo.) a 


-0.07 ± 0.07 


0.70 ± 0.10 


0.21 ± 0.08 




this work (no evo.) 


-0.11 ± 0.08 


0.67 ± 0.10 


0.23 ± 0.09 




Bcntz ct al. (with evo.) 


0.06 ± 0.06 


0.72 ± 0.09 


0.20 ± 0.06 




Bcntz ct al. (no evo.) 


0.02 ± 0.06 


0.70 ± 0.08 


0.20 ± 0.06 


BCES 


this work (with evo.) 


-0.12 ± 0.06 


0.81 ± 0.11 






this work (no evo.) 


-0.15 ± 0.06 


0.77 ± 0.10 






Bcntz ct al. (with evo.) 


0.02 ± 0.06 


0.84 ± 0.09 






Bentz et al. (no evo.) 


-0.02 ± 0.06 


0.80 ± 0.09 




linear fit 


this work (host; with cvo.) a 


-0.38 ± 0.12 


0.96 ± 0.18 


0.24 ± 0.11 



NOTE. — Comparison between the different fits (in the form of Equation 5) to the local 
RM AGN log A/bh - log £ B ph,v relation: with and without correction for passive lumi- 
nosity evolution, different fitting methods (linear fit vs BCES), and Bcntz ct al. (2009b) 
results vs. new analysis in this paper. In the last row, we give the fit to the local RM 
AGN log A/bh - log £host,v relation derived in this paper. 

Col. (1): Fitting method. Linear fit with intrinsic scatter or BCES for comparison with 
Bcntz ct al. (2009a). Col. (2): Sample, "evo." indicates whether or not data have 
been corrected for luminosity evolution. Col. (3): Mean and uncertainty on the best fit 
intercept. Col. (4): Mean and uncertainty on the best fit slope. Col. (5): Mean and 
uncertainty on the best fit intrinsic scatter (for "linear fit" only). 
a This is the fit we use in the subsequent analysis. 
b This is the fit used in Bcn tz et all d2009aD . 
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